home *** CD-ROM | disk | FTP | other *** search
/ The Atari Compendium / The Atari Compendium (Toad Computers) (1994).iso / files / prgtools / gnustuff / tos / updates / update31.zoo / lib / diffs < prev   
Encoding:
Text File  |  1993-05-29  |  90.5 KB  |  3,663 lines
  1. ===================================================================
  2. RCS file: /net/acae127/home/bammi/etc/src/master/atari/lib/Changelo,v
  3. retrieving revision 1.93
  4. diff -c -r1.93 Changelo
  5. *** 1.93    1993/05/21 12:52:46
  6. --- Changelo    1993/05/29 20:26:25
  7. ***************
  8. *** 3730,3732 ****
  9. --- 3730,3767 ----
  10.       Fixed sign for retinf, and retzero.
  11.   
  12.   ---------------------------- Patchlevel 89 ---------------------------
  13. + 
  14. +     michal did a tremendous amount of job tracking down modf and
  15. +     fixing the following. thanks michal.
  16. + 
  17. +     Wed May 26 20:35:34 1993  Michal Jaegermann  (michal at smok)
  18. + 
  19. +     * fixed norm_df and norm_sf not to set errno when processing
  20. +       legitimate zero argument
  21. +     * fixed wrong error codes and made sure that future changes
  22. +       in errno.h will be automatically reflected;
  23. +       this involves changes in errbase.h and errno.h
  24. +     * removed unused code from _divsf3.cpp
  25. +     * corrected bug in setting sign of Inf in _truncdf.cpp
  26. +     * made zeros returned from frexp and ldexp consistenly signed
  27. +       (they were signed or unsigned, depending on a way they were
  28. +        calculated)
  29. +     * renamed _cmpdf2.s to _cmpdf2.cpp and _cmpsf2.s to _cmpsf2.cpp
  30. +       and added preprocessor directives to split them into separate
  31. +       objects while compiling (required change in mincl)
  32. +     * replaced 68000 code for modf() with totally new, shorter and
  33. +       much faster version
  34. +     * added Olaf Flebbe fixes for sign of 0 in multiplication routines
  35. +     * general code cleanup all over the place in floating point support
  36. +       routines
  37. + 
  38. + scanf.c:: ++jrb
  39. +     fix bug reported by warwick, when "%Ns", the test for the end of
  40. +     the %s loop was incorrect. (left off 1 char too soon).
  41. + 
  42. +     added some decls for non ansi platforms.
  43. + 
  44. +     use long double type only when __M68881__ not with soft float.
  45. +     out soft float stuff knows nothing about long doubles.
  46. + 
  47. + ---------------------------- Patchlevel 90 ---------------------------
  48. ===================================================================
  49. RCS file: /net/acae127/home/bammi/etc/src/master/atari/lib/PatchLev.h,v
  50. retrieving revision 1.65
  51. diff -c -r1.65 PatchLev.h
  52. *** 1.65    1993/05/21 12:52:51
  53. --- PatchLev.h    1993/05/29 20:26:29
  54. ***************
  55. *** 1,5 ****
  56.   
  57. ! #define    PatchLevel "89"
  58.   
  59.   /*
  60.    *
  61. --- 1,5 ----
  62.   
  63. ! #define    PatchLevel "90"
  64.   
  65.   /*
  66.    *
  67. ===================================================================
  68. RCS file: /net/acae127/home/bammi/etc/src/master/atari/lib/fscanf.c,v
  69. retrieving revision 1.9
  70. diff -c -r1.9 fscanf.c
  71. *** 1.9    1992/01/29 20:58:29
  72. --- fscanf.c    1993/05/29 20:26:32
  73. ***************
  74. *** 5,11 ****
  75.   #include <compiler.h>
  76.   #endif
  77.   
  78. ! #ifdef sun
  79.   # define fungetc ungetc
  80.   #endif
  81.   
  82. --- 5,11 ----
  83.   #include <compiler.h>
  84.   #endif
  85.   
  86. ! #if defined(sun) || defined(__hpux)
  87.   # define fungetc ungetc
  88.   #endif
  89.   
  90. ===================================================================
  91. RCS file: /net/acae127/home/bammi/etc/src/master/atari/lib/mincl,v
  92. retrieving revision 1.27
  93. diff -c -r1.27 mincl
  94. *** 1.27    1993/05/21 12:52:54
  95. --- mincl    1993/05/29 20:26:34
  96. ***************
  97. *** 1,6 ****
  98.   #
  99.   # GCC specific stuff
  100. ! GLIB1 = _addsubd.o _addsubs.o _cmpdf2.o _cmpsf2.o _divdf3.o _divmods.o \
  101.       _divsf3.o _extends.o _fixdfsi.o _fxunsd.o _fltsi.o \
  102.       _muldf3.o _mulsf3.o _mulsi3.o _negdf2.o _negsf2.o _normdf.o _normsf.o \
  103.       _truncdf.o _udivmod.o _umulsi3.o _fixsfsi.o _fltsisf.o _isnan.o
  104. --- 1,9 ----
  105.   #
  106.   # GCC specific stuff
  107. ! FCOMPD = _eqdf2.o _gtdf2.o _ledf2.o _gedf2.o _ltdf2.o
  108. ! FCOMPS = _eqsf2.o _gtsf2.o _lesf2.o _gesf2.o _ltsf2.o
  109. ! 
  110. ! GLIB1 = _addsubd.o _addsubs.o _divdf3.o _divmods.o \
  111.       _divsf3.o _extends.o _fixdfsi.o _fxunsd.o _fltsi.o \
  112.       _muldf3.o _mulsf3.o _mulsi3.o _negdf2.o _negsf2.o _normdf.o _normsf.o \
  113.       _truncdf.o _udivmod.o _umulsi3.o _fixsfsi.o _fltsisf.o _isnan.o
  114. ***************
  115. *** 12,18 ****
  116.       _caps_New.o _builtin_del.o _trampoline.o __main.o _ctor_list.o \
  117.       _dtor_list.o _ffsdi2.o
  118.   
  119. ! GCC=     $(GLIB1) $(GLIB2) \
  120.       ldexp.o frexp.o modf.o alloca.o setjmp.o osbind.o\
  121.       linea.o  alglobal.o sysvar.o gmon.o screen.o stksiz.o binmode.o \
  122.       bblink.o defmode.o
  123. --- 15,21 ----
  124.       _caps_New.o _builtin_del.o _trampoline.o __main.o _ctor_list.o \
  125.       _dtor_list.o _ffsdi2.o
  126.   
  127. ! GCC=     $(GLIB1) $(GLIB2) $(FCOMPD) $(FCOMPS)\
  128.       ldexp.o frexp.o modf.o alloca.o setjmp.o osbind.o\
  129.       linea.o  alglobal.o sysvar.o gmon.o screen.o stksiz.o binmode.o \
  130.       bblink.o defmode.o
  131. ***************
  132. *** 86,91 ****
  133. --- 89,100 ----
  134.   #
  135.   # Assembler stuff that needs to be preprocessed by cpp
  136.   #
  137. + $(FCOMPD): %.o: _cmpdf2.cpp
  138. +     $(CC) $(PPFLAGS) -DL$* -c $< -o $@
  139. + 
  140. + $(FCOMPS): %.o: _cmpsf2.cpp
  141. +     $(CC) $(PPFLAGS) -DL$* -c $< -o $@
  142. + 
  143.   %.o : %.cpp
  144.       $(CC) $(PPFLAGS) -c $<
  145.   
  146. ===================================================================
  147. RCS file: /net/acae127/home/bammi/etc/src/master/atari/lib/scanf.c,v
  148. retrieving revision 1.15
  149. diff -c -r1.15 scanf.c
  150. *** 1.15    1992/10/09 20:32:44
  151. --- scanf.c    1993/05/29 20:26:36
  152. ***************
  153. *** 39,44 ****
  154. --- 39,45 ----
  155.   #include <stdio.h>
  156.   #include <stdlib.h>
  157.   #include <string.h>
  158. + #include <compiler.h>
  159.   
  160.   /* the code assumes this definition, note: traditional <ctype> def
  161.    * of tolower will break the code. Ansi def should be ok.
  162. ***************
  163. *** 70,75 ****
  164. --- 71,78 ----
  165.   #define input_error()    return( done < 1 ? EOF : done )
  166.   #define    memory_error()    return((errno = ENOMEM), EOF)
  167.   
  168. + __EXTERN long int strtol __PROTO((const char *nptr, char **endptr, int base));
  169. + __EXTERN double strtod __PROTO((const char *s, char **endptr));
  170.   
  171.   /* Read formatted input from S according to the format string
  172.      FORMAT, using the argument list in ARG.
  173. ***************
  174. *** 110,116 ****
  175. --- 113,123 ----
  176.       unsigned long int unum;
  177.   #if FLOATS
  178.       /* Floating-point holding variable.  */
  179. + # ifdef __M68881__
  180.       long double fp_num;
  181. + # else
  182. +     double fp_num;
  183. + # endif
  184.   #endif
  185.       /* Character-buffer pointer.  */
  186.       register char *str;
  187. ***************
  188. *** 285,293 ****
  189.               break;
  190.           if (do_assign)
  191.               *str++ = c;
  192. !         if (width > 0 && --width == 0)
  193. !             break;
  194. !         } while (inchar() != EOF);
  195.           
  196.           if (do_assign)
  197.           {
  198. --- 292,298 ----
  199.               break;
  200.           if (do_assign)
  201.               *str++ = c;
  202. !         } while ((inchar() != EOF) && ((width > 0) ? --width != 0 : 1));
  203.           
  204.           if (do_assign)
  205.           {
  206. ***************
  207. *** 468,474 ****
  208. --- 473,483 ----
  209.           if (do_assign)
  210.           {
  211.           if (is_long_double)
  212. + #ifdef __M68881__
  213.               *va_arg(arg, long double *) = fp_num;
  214. + #else
  215. +             *va_arg(arg, double *) = fp_num;
  216. + #endif
  217.           else if (is_long)
  218.               *va_arg(arg, double *) = (double) fp_num;
  219.           else
  220. ===================================================================
  221. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_addsubd.cpp,v
  222. retrieving revision 1.6
  223. diff -c -r1.6 _addsubd.cpp
  224. *** 1.6    1993/03/01 17:33:55
  225. --- _addsubd.cpp    1993/05/29 20:28:20
  226. ***************
  227. *** 10,15 ****
  228. --- 10,19 ----
  229.   |
  230.   | written by Kai-Uwe Bloem (I5110401@dbstu1.bitnet).
  231.   | Based on a 80x86 floating point packet from comp.os.minix, written by P.Housel
  232. + |
  233. + | Revision 1.3.5 michal 05-93 (ntomczak@vm.ucs.ualberta.ca)
  234. + |   + code smoothing
  235. + |
  236.   | patched by Olaf Flebbe (flebbe@tat.physik.uni-tuebingen.de)
  237.   |
  238.   | Revision 1.3.4 olaf 11-92 :
  239. ***************
  240. *** 23,30 ****
  241.   |
  242.   | Revision 1.3.3 olaf 11-92 :
  243.   |  + changed to get rid of rounding bits. a sticky register (d3) is
  244. ! |    sufficient. 
  245. ! | 
  246.   | Revision 1.3.2 olaf 10-92 :
  247.   |  + increased comparson by one again. (Dont understand, but it works)
  248.   |  + corrected negation of rounding bits and mantissa
  249. --- 27,34 ----
  250.   |
  251.   | Revision 1.3.3 olaf 11-92 :
  252.   |  + changed to get rid of rounding bits. a sticky register (d3) is
  253. ! |    sufficient.
  254. ! |
  255.   | Revision 1.3.2 olaf 10-92 :
  256.   |  + increased comparson by one again. (Dont understand, but it works)
  257.   |  + corrected negation of rounding bits and mantissa
  258. ***************
  259. *** 59,85 ****
  260.       moveml    d2-d7,sp@-    | save registers
  261.       moveml    a0@,d4-d5/d6-d7    | d4-d5 = v, d6-d7 = u
  262.   
  263.       movel    d6,d0        | d0 = u.exp
  264.       swap    d0
  265. !     movel    d6,d2        | d2.h = u.sign
  266. !     movew    d0,d2
  267.       lsrw    #4,d0
  268. -     andw    #0x07ff,d0    | kill sign bit
  269.   
  270.       movel    d4,d1        | d1 = v.exp
  271.       swap    d1
  272. !     eorw    d1,d2        | d2.l = u.sign ^ v.sign
  273.       lsrw    #4,d1
  274. -     andw    #0x07ff,d1    | kill sign bit
  275. - 
  276. -     andl    #0x0fffff,d6    | remove exponent from u.mantissa
  277. -     andl    #0x0fffff,d4    | remove exponent from v.mantissa
  278.   |
  279.   | Now perform testing of NaN And infinities
  280.   |
  281. !     cmpw    #0x07ff,d0
  282.       beq    0f
  283. !     cmpw    #0x07ff,d1
  284.       bne    nospec
  285.       bra    1f
  286.   |    first operand is special
  287. --- 63,90 ----
  288.       moveml    d2-d7,sp@-    | save registers
  289.       moveml    a0@,d4-d5/d6-d7    | d4-d5 = v, d6-d7 = u
  290.   
  291. +     movel    #0x0fffff,d3
  292.       movel    d6,d0        | d0 = u.exp
  293. +     andl    d3,d6        | remove exponent from u.mantissa
  294. +     movel    d0,d2        | d2.h = u.sign
  295.       swap    d0
  296. !     movew    d0,d2        | d2.l = u.sign
  297. !     bclr    #15,d0        | kill sign bit
  298.       lsrw    #4,d0
  299.   
  300.       movel    d4,d1        | d1 = v.exp
  301. +     andl    d3,d4        | remove exponent from v.mantissa
  302.       swap    d1
  303. !     eorw    d1,d2        | d2.l = u.sign ^ v.sign (in bit 15)
  304. !     bclr    #15,d1        | kill sign bit
  305.       lsrw    #4,d1
  306.   |
  307.   | Now perform testing of NaN And infinities
  308.   |
  309. !     movew    #0x7ff,d3
  310. !     cmpw    d3,d0
  311.       beq    0f
  312. !     cmpw    d3,d1
  313.       bne    nospec
  314.       bra    1f
  315.   |    first operand is special
  316. ***************
  317. *** 88,100 ****
  318.       beq    bothspec
  319.   |
  320.   |    u is special
  321. ! | 
  322.       movel    d7,d0
  323.       orl    d6,d0
  324.       bne    retnan        | arith with Nan gives always Nan
  325.   
  326.       movel    a0@(8),d0        | copy infinity
  327. !     clrl    d1
  328.       bra     return
  329.   |
  330.   |    v is special
  331. --- 93,105 ----
  332.       beq    bothspec
  333.   |
  334.   |    u is special
  335. ! |
  336.       movel    d7,d0
  337.       orl    d6,d0
  338.       bne    retnan        | arith with Nan gives always Nan
  339.   
  340.       movel    a0@(8),d0        | copy infinity
  341. !     moveq    #0,d1
  342.       bra     return
  343.   |
  344.   |    v is special
  345. ***************
  346. *** 103,109 ****
  347.       orl    d4,d0
  348.       bne    retnan
  349.       movel    a0@,d0
  350. !     clrl    d1
  351.       bra    return
  352.   |
  353.   |    u and v are both special
  354. --- 108,114 ----
  355.       orl    d4,d0
  356.       bne    retnan
  357.       movel    a0@,d0
  358. !     moveq    #0,d1
  359.       bra    return
  360.   |
  361.   |    u and v are both special
  362. ***************
  363. *** 113,156 ****
  364.       orl    d6,d0
  365.       orl    d5,d0
  366.       orl    d4,d0
  367. !     beq    bothinf
  368. !     bra    retnan    
  369.   |
  370.   | Both are infinities Test if cancellation
  371. ! | 
  372.   bothinf:
  373.       tstw    d2
  374. !     bpl    retinf    
  375.   |
  376.   | return a quiet NaN
  377.   |
  378. ! retnan:    movel    #0x7fffffff,d0
  379. !     moveql  #-1,d1
  380.       bra    return
  381.   
  382. ! retinf: movel   #0x7ff00000,d0
  383. !     clrl    d1
  384. !     tstl    d2
  385. !     bpl    return
  386. !     bchg    #31,d0
  387.   return:    moveml    sp@+,d2-d7
  388.       rts
  389. ! |    
  390.   | Ok, no inifinty or Nan involved..
  391.   |
  392. ! nospec:    tstw    d0        | check for zero exponent - no leading "1"
  393. !     beq    0f
  394. !     bset    #20,d6        | restore implied leading "1"
  395. !     bra    1f
  396. ! 0:    addw    #1,d0        | "normalize" exponent
  397.   1:
  398. -     tstw    d1        | check for zero exponent - no leading "1"
  399. -     beq    0f
  400.       bset    #20,d4        | restore implied leading "1"
  401. !     bra    1f
  402. ! 0:    addw    #1,d1        | "normalize" exponent
  403.   1:
  404. !     clrl    d3        | init sticky register
  405.       negw    d1        | d1 = u.exp - v.exp
  406.       addw    d0,d1
  407.       beq    5f        | exponents are equal - no shifting neccessary
  408. --- 118,160 ----
  409.       orl    d6,d0
  410.       orl    d5,d0
  411.       orl    d4,d0
  412. !     bne    retnan
  413.   |
  414.   | Both are infinities Test if cancellation
  415. ! |
  416.   bothinf:
  417.       tstw    d2
  418. !     bpl    retinf
  419.   |
  420.   | return a quiet NaN
  421.   |
  422. ! retnan: moveql    #-1,d1
  423. !     movel    d1,d0
  424. !     bclr    #31,d0        | 0x7fffffff -> d0
  425.       bra    return
  426.   
  427. ! retinf:    moveq    #0,d1
  428. !     movel    #0xffe00000,d0    | we will return #0xfff00000 or #0x7ff00000
  429. !     lslw    #1,d2
  430. !     roxrl   #1,d0        | shift in high bit as given by d2
  431.   return:    moveml    sp@+,d2-d7
  432.       rts
  433. ! |
  434.   | Ok, no inifinty or Nan involved..
  435.   |
  436. ! nospec:    bset    #20,d6        | restore implied leading "1"
  437. !     tstw    d0        | check for zero exponent - no leading "1"
  438. !     bne    1f
  439. !     bclr    #20,d6        | no implied leading "1", instead ...
  440. !     addw    #1,d0        | "normalize" exponent
  441.   1:
  442.       bset    #20,d4        | restore implied leading "1"
  443. !     tstw    d1        | check for zero exponent - no leading "1"
  444. !     bne    1f
  445. !     bclr    #20,d4        | no implied leading "1", instead ...
  446. !     addw    #1,d1        | "normalize" exponent
  447.   1:
  448. !     moveq    #0,d3        | init sticky register
  449.       negw    d1        | d1 = u.exp - v.exp
  450.       addw    d0,d1
  451.       beq    5f        | exponents are equal - no shifting neccessary
  452. ***************
  453. *** 170,186 ****
  454.   | most significant digit, while gaining an additional digit for
  455.   | rounding.
  456.   |
  457. !     moveql  #1,d3    
  458.   2:    addl    d7,d7
  459.       addxl    d6,d6
  460.       subw    #1,d0        | decrement exponent
  461.       subw    #1,d1        | decrement counter
  462.       dbeq    d3,2b
  463. !     clrl    d3
  464.   |
  465.   | now shift other mantissa right as fast as possible (almost).
  466.   |
  467. ! 3:    
  468.       cmpw    #16,d1        | see if fast rotate possible
  469.       blt    4f
  470.       orw    d5,d3        | set sticky word
  471. --- 174,190 ----
  472.   | most significant digit, while gaining an additional digit for
  473.   | rounding.
  474.   |
  475. !     moveql  #1,d3
  476.   2:    addl    d7,d7
  477.       addxl    d6,d6
  478.       subw    #1,d0        | decrement exponent
  479.       subw    #1,d1        | decrement counter
  480.       dbeq    d3,2b
  481. !     moveq    #0,d3
  482.   |
  483.   | now shift other mantissa right as fast as possible (almost).
  484.   |
  485. ! 3:
  486.       cmpw    #16,d1        | see if fast rotate possible
  487.       blt    4f
  488.       orw    d5,d3        | set sticky word
  489. ***************
  490. *** 193,203 ****
  491.   
  492.   0:    moveb   d5,d2        | use d2.b as scratch
  493.       andb    #1,d2        | test if 1 is shifted out
  494. !     orb    d2,d3        | and put it in sticky 
  495.       lsrl    #1,d4        | shift v.mant right the rest of the way
  496.       roxrl    #1,d5        | to line it up with u.mant
  497.   4:    dbra    d1,0b        | loop
  498. !     
  499.   5:
  500.       tstw    d2        | are the signs equal ?
  501.       bpl    6f        | yes, no negate necessary
  502. --- 197,207 ----
  503.   
  504.   0:    moveb   d5,d2        | use d2.b as scratch
  505.       andb    #1,d2        | test if 1 is shifted out
  506. !     orb    d2,d3        | and put it in sticky
  507.       lsrl    #1,d4        | shift v.mant right the rest of the way
  508.       roxrl    #1,d5        | to line it up with u.mant
  509.   4:    dbra    d1,0b        | loop
  510. ! 
  511.   5:
  512.       tstw    d2        | are the signs equal ?
  513.       bpl    6f        | yes, no negate necessary
  514. ***************
  515. *** 205,218 ****
  516.   | negate secand mantissa. One has to check the sticky word in order
  517.   | to correct the twos complement.
  518.   |
  519. !     tstw    d3        | 
  520.       beq     9f        | No cerrection necessary
  521. !     clrl    d1
  522.       addql   #1,d5
  523.       addxl   d1,d4
  524.   9:    negl    d5
  525.       negxl    d4
  526. !         
  527.   6:
  528.       addl    d5,d7        | u.mant = u.mant + v.mant
  529.       addxl    d4,d6
  530. --- 209,222 ----
  531.   | negate secand mantissa. One has to check the sticky word in order
  532.   | to correct the twos complement.
  533.   |
  534. !     tstw    d3        |
  535.       beq     9f        | No cerrection necessary
  536. !     moveq    #0,d1
  537.       addql   #1,d5
  538.       addxl   d1,d4
  539.   9:    negl    d5
  540.       negxl    d4
  541. ! 
  542.   6:
  543.       addl    d5,d7        | u.mant = u.mant + v.mant
  544.       addxl    d4,d6
  545. ***************
  546. *** 226,232 ****
  547.   7:
  548.       movel    d6,d4        | move result for normalization
  549.       movel    d7,d5
  550. !     clrl    d1
  551.       tstl    d3
  552.       beq     8f
  553.       moveql   #-1,d1
  554. --- 230,236 ----
  555.   7:
  556.       movel    d6,d4        | move result for normalization
  557.       movel    d7,d5
  558. !     moveq    #0,d1
  559.       tstl    d3
  560.       beq     8f
  561.       moveql   #-1,d1
  562. ***************
  563. *** 267,278 ****
  564.       lea    0xfffffa50:w,a0
  565.       movew    #0x5400,a0@(comm)    | load first argument to fp0
  566.       cmpiw    #0x8900,a0@(resp)    | check
  567. !     movel    a7@(4),a0@
  568. !     movel    a7@(8),a0@
  569.       movew    #0x5428,a0@(comm)
  570.       .long    0x0c688900, 0xfff067f8
  571. !     movel    a7@(12),a0@
  572. !     movel    a7@(16),a0@
  573.       movew    #0x7400,a0@(comm)    | result to d0/d1
  574.       .long    0x0c688900, 0xfff067f8
  575.       movel    a0@,d0
  576. --- 271,282 ----
  577.       lea    0xfffffa50:w,a0
  578.       movew    #0x5400,a0@(comm)    | load first argument to fp0
  579.       cmpiw    #0x8900,a0@(resp)    | check
  580. !     movel    sp@(4),a0@
  581. !     movel    sp@(8),a0@
  582.       movew    #0x5428,a0@(comm)
  583.       .long    0x0c688900, 0xfff067f8
  584. !     movel    sp@(12),a0@
  585. !     movel    sp@(16),a0@
  586.       movew    #0x7400,a0@(comm)    | result to d0/d1
  587.       .long    0x0c688900, 0xfff067f8
  588.       movel    a0@,d0
  589. ***************
  590. *** 284,295 ****
  591.       lea    0xfffffa50:w,a0
  592.       movew    #0x5400,a0@(comm)        | load fp0
  593.       cmpiw    #0x8900,a0@(resp)        | got it?
  594. !     movel    a7@(4),a0@            | take a hi from stack to FPU
  595. !     movel    a7@(8),a0@            | take a lo from stack to FPU
  596.       movew    #0x5422,a0@(comm)        | add second arg to fp0
  597.       .long    0x0c688900, 0xfff067f8
  598. !     movel    a7@(12),a0@            | move b hi from stack to FPU
  599. !     movel    a7@(16),a0@            | move b lo from stack to FPU
  600.       movew    #0x7400,a0@(comm)        | result to d0/d1
  601.       .long    0x0c688900, 0xfff067f8
  602.       movel    a0@,d0                | download result
  603. --- 288,299 ----
  604.       lea    0xfffffa50:w,a0
  605.       movew    #0x5400,a0@(comm)        | load fp0
  606.       cmpiw    #0x8900,a0@(resp)        | got it?
  607. !     movel    sp@(4),a0@            | take a hi from stack to FPU
  608. !     movel    sp@(8),a0@            | take a lo from stack to FPU
  609.       movew    #0x5422,a0@(comm)        | add second arg to fp0
  610.       .long    0x0c688900, 0xfff067f8
  611. !     movel    sp@(12),a0@            | move b hi from stack to FPU
  612. !     movel    sp@(16),a0@            | move b lo from stack to FPU
  613.       movew    #0x7400,a0@(comm)        | result to d0/d1
  614.       .long    0x0c688900, 0xfff067f8
  615.       movel    a0@,d0                | download result
  616. ===================================================================
  617. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_addsubs.cpp,v
  618. retrieving revision 1.6
  619. diff -c -r1.6 _addsubs.cpp
  620. *** 1.6    1993/03/01 17:33:57
  621. --- _addsubs.cpp    1993/05/29 20:28:21
  622. ***************
  623. *** 10,15 ****
  624. --- 10,19 ----
  625.   |
  626.   | written by Kai-Uwe Bloem (I5110401@dbstu1.bitnet).
  627.   | Based on a 80x86 floating point packet from comp.os.minix, written by P.Housel
  628. + |
  629. + | Revision 1.3.5 michal 05-93 (ntomczak@vm.ucs.ualberta.ca)
  630. + |   + code smoothing
  631. + |
  632.   | patched by Olaf Flebbe (flebbe@tat.physik.uni-tuebingen.de)
  633.   |
  634.   | Revision 1.3.4 olaf 11-92 :
  635. ***************
  636. *** 61,88 ****
  637.       moveml    d2-d5,sp@-    | save registers
  638.       moveml    a0@,d4/d5    | d4 = v, d5 = u
  639.   
  640.       movel    d5,d0        | d0 = u.exp
  641. -     swap    d0
  642.       movel    d5,d2        | d2.h = u.sign
  643. !     movew    d0,d2
  644.       lsrw    #7,d0
  645. -     andw    #0xff,d0    | kill sign bit (exponent is 8 bits)
  646.   
  647.       movel    d4,d1        | d1 = v.exp
  648.       swap    d1
  649. !     eorw    d1,d2        | d2.l = u.sign ^ v.sign
  650.       lsrw    #7,d1
  651. -     andw    #0xff,d1    | kill sign bit (exponent is 8 bits)
  652. - 
  653. -     andl    #0x7fffff,d5    | remove exponent from mantissa
  654. -     andl    #0x7fffff,d4    | remove exponent from mantissa
  655.   |
  656.   | Now perform testing of NaN And infinities
  657.   |
  658. !     cmpb    #0xff,d0
  659.       beq    0f
  660. !     cmpb    #0xff,d1
  661. !     bne    nospec        
  662.       bra    1f
  663.   |    first operand is special
  664.   |
  665. --- 65,93 ----
  666.       moveml    d2-d5,sp@-    | save registers
  667.       moveml    a0@,d4/d5    | d4 = v, d5 = u
  668.   
  669. +     movel    #0x7fffff,d3
  670.       movel    d5,d0        | d0 = u.exp
  671.       movel    d5,d2        | d2.h = u.sign
  672. !     swap    d0
  673. !     movew    d0,d2        | d2 = u.sign
  674. !     andl    d3,d5        | remove exponent from u.mantissa
  675. !     bclr    #15,d0        | kill sign bit
  676.       lsrw    #7,d0
  677.   
  678.       movel    d4,d1        | d1 = v.exp
  679. +     andl    d3,d4        | remove exponent from v.mantissa
  680.       swap    d1
  681. !     eorw    d1,d2        | d2 = u.sign ^ v.sign (in bit 15)
  682. !     bclr    #15,d1        | kill sign bit
  683.       lsrw    #7,d1
  684.   |
  685.   | Now perform testing of NaN And infinities
  686.   |
  687. !     moveq    #-1,d3
  688. !     cmpb    d3,d0
  689.       beq    0f
  690. !     cmpb    d3,d1
  691. !     bne    nospec
  692.       bra    1f
  693.   |    first operand is special
  694.   |
  695. ***************
  696. *** 90,100 ****
  697.       beq    bothspec
  698.   |
  699.   |    u is special
  700. ! | 
  701.       tstl    d5
  702.       bne    retnan        | arith with Nan gives always Nan
  703.   
  704. !     movel    a0@(4),d0        | copy infinity
  705.       bra     return
  706.   |
  707.   |    v is special
  708. --- 95,105 ----
  709.       beq    bothspec
  710.   |
  711.   |    u is special
  712. ! |
  713.       tstl    d5
  714.       bne    retnan        | arith with Nan gives always Nan
  715.   
  716. !     movel    a0@(4),d0    | copy infinity with its sign
  717.       bra     return
  718.   |
  719.   |    v is special
  720. ***************
  721. *** 109,150 ****
  722.   bothspec:
  723.       movel    d5,d0
  724.       orl    d4,d0
  725. !     beq    bothinf
  726. !     bra    retnan    
  727.   |
  728.   | Both are infinities Test if cancellation
  729. ! | 
  730.   bothinf:
  731.       tstw    d2
  732. !     bpl    retinf    
  733.   |
  734.   | return a quiet NaN
  735.   |
  736. ! retnan:    movel    #0x7fffffff,d0
  737.       bra    return
  738.   
  739. ! retinf: movel   #0x7f800000,d0
  740. !     tstl    d2
  741. !     bpl    return
  742. !     bchg    #31,d0
  743.   return:    moveml    sp@+,d2-d5
  744.       rts
  745. ! |    
  746.   | Ok, no inifinty or Nan involved..
  747.   |
  748. ! nospec:    tstw    d0        | check for zero exponent - no leading "1"
  749. !     beq    0f
  750. !     bset    #23,d5        | restore implied leading "1"
  751. !     bra    1f
  752. ! 0:    addw    #1,d0        | "normalize" exponent
  753.   1:
  754. -     tstw    d1        | check for zero exponent - no leading "1"
  755. -     beq    0f
  756.       bset    #23,d4        | restore implied leading "1"
  757. !     bra    1f
  758. ! 0:    addw    #1,d1        | "normalize" exponent
  759.   1:
  760. !     clrl    d3        | (put initial zero rounding bits in d3)
  761.       negw    d1        | d1 = u.exp - v.exp
  762.       addw    d0,d1
  763.       beq    5f        | exponents are equal - no shifting neccessary
  764. --- 114,154 ----
  765.   bothspec:
  766.       movel    d5,d0
  767.       orl    d4,d0
  768. !     bne    retnan
  769.   |
  770.   | Both are infinities Test if cancellation
  771. ! |
  772.   bothinf:
  773.       tstw    d2
  774. !     bpl    retinf
  775.   |
  776.   | return a quiet NaN
  777.   |
  778. ! retnan: moveql    #-1,d0
  779. !     lsrl    #1,d0        | 0x7fffffff -> d0
  780.       bra    return
  781.   
  782. ! retinf: movel    #0xff000000,d0  | we will return #0xff800000 or #0x7f800000
  783. !     lslw    #1,d2
  784. !     roxrl   #1,d0        | shift in high bit as given by d2
  785.   return:    moveml    sp@+,d2-d5
  786.       rts
  787. ! |
  788.   | Ok, no inifinty or Nan involved..
  789.   |
  790. ! nospec:    bset    #23,d5        | restore implied leading "1"
  791. !     tstw    d0        | check for zero exponent - no leading "1"
  792. !     bne    1f
  793. !     bclr    #23,d5        | remove it
  794. !     addqw    #1,d0        | "normalize" exponent
  795.   1:
  796.       bset    #23,d4        | restore implied leading "1"
  797. !     tstw    d1        | check for zero exponent - no leading "1"
  798. !     bne    1f
  799. !     bclr    #23,d4        | remove it
  800. !     addqw    #1,d1        | "normalize" exponent
  801.   1:
  802. !     moveq    #0,d3        | (put initial zero rounding bits in d3)
  803.       negw    d1        | d1 = u.exp - v.exp
  804.       addw    d0,d1
  805.       beq    5f        | exponents are equal - no shifting neccessary
  806. ***************
  807. *** 163,174 ****
  808.   | most significant digit, while gaining an additional digit for
  809.   | rounding.
  810.   |
  811. !     moveql    #1,d3    
  812.   2:    addl    d5,d5
  813. !     subw    #1,d0        | decrement exponent
  814. !     subw    #1,d1        | done shifting altogether ?
  815.       dbeq    d3,2b        | loop if still can shift u.mant more
  816. !     clrl    d3
  817.   
  818.       cmpw    #16,d1        | see if fast rotate possible
  819.       blt    4f
  820. --- 167,178 ----
  821.   | most significant digit, while gaining an additional digit for
  822.   | rounding.
  823.   |
  824. !     moveql    #1,d3
  825.   2:    addl    d5,d5
  826. !     subqw    #1,d0        | decrement exponent
  827. !     subqw    #1,d1        | done shifting altogether ?
  828.       dbeq    d3,2b        | loop if still can shift u.mant more
  829. !     moveq    #0,d3
  830.   
  831.       cmpw    #16,d1        | see if fast rotate possible
  832.       blt    4f
  833. ===================================================================
  834. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_divdf3.cpp,v
  835. retrieving revision 1.6
  836. diff -c -r1.6 _divdf3.cpp
  837. *** 1.6    1993/05/21 12:55:07
  838. --- _divdf3.cpp    1993/05/29 20:28:23
  839. ***************
  840. *** 23,34 ****
  841.   
  842.       pea    pc@(LC0)
  843.       pea    Stderr
  844. !     jbsr    _fprintf    |
  845. !     addql    #8,a7        |
  846. !                     | set _errno to ERANGE
  847. !     moveq    #ERANGE,d0
  848.       Emove    d0,Errno
  849. !     moveml    __infinitydf,d0-d1    | return signed infinity 
  850.       btst    #31,a7@(4)        | transfer sign of dividend
  851.       beq    clear            | (mjr++)
  852.       bset    #31,d0            |
  853. --- 23,34 ----
  854.   
  855.       pea    pc@(LC0)
  856.       pea    Stderr
  857. !     jbsr    _fprintf
  858. !     addql    #8,a7
  859. ! 
  860. !     moveq    #Erange,d0        | set _errno to ERANGE
  861.       Emove    d0,Errno
  862. !     moveml    __infinitydf,d0-d1    | return signed infinity
  863.       btst    #31,a7@(4)        | transfer sign of dividend
  864.       beq    clear            | (mjr++)
  865.       bset    #31,d0            |
  866. ***************
  867. *** 44,53 ****
  868.   | written by Kai-Uwe Bloem (I5110401@dbstu1.bitnet).
  869.   | Based on a 80x86 floating point packet from comp.os.minix, written by P.Housel
  870.   |
  871.   | patched by Olaf Flebbe (flebbe@tat.physik.uni-tuebingen.de)
  872.   |
  873.   | Revision 1.2.3 olaf 4-93
  874. ! |   + Fixed sign for retinf, and retzero: it is in d2.w 
  875.   |
  876.   | Revision 1.2.2 olaf 12-92
  877.   |   + added support for NaN and Infinites
  878. --- 44,57 ----
  879.   | written by Kai-Uwe Bloem (I5110401@dbstu1.bitnet).
  880.   | Based on a 80x86 floating point packet from comp.os.minix, written by P.Housel
  881.   |
  882. + | Revision 1.2.4 michal 05-93 (ntomczak@vm.ucs.ualberta.ca)
  883. + |   + resynchro with errno codes
  884. + |   + code smoothing
  885. + |
  886.   | patched by Olaf Flebbe (flebbe@tat.physik.uni-tuebingen.de)
  887.   |
  888.   | Revision 1.2.3 olaf 4-93
  889. ! |   + Fixed sign for retinf, and retzero: it is in d2.w
  890.   |
  891.   | Revision 1.2.2 olaf 12-92
  892.   |   + added support for NaN and Infinites
  893. ***************
  894. *** 73,98 ****
  895.       moveml    d2-d7,sp@-    | save registers
  896.       moveml    a0@,d4-d5/d6-d7    | d4-d5 = u, d6-d7 = v
  897.   
  898.       movel    d4,d0        | d0 = u.exp
  899.       swap    d0
  900.       movew    d0,d2        | d2 = u.sign
  901.       lsrw    #4,d0
  902. -     andw    #0x07ff,d0    | kill sign bit
  903.   
  904.       movel    d6,d1        | d1 = v.exp
  905.       swap    d1
  906.       eorw    d1,d2        | d2 = u.sign ^ v.sign (in bit 15)
  907.       lsrw    #4,d1
  908. -     andw    #0x07ff,d1    | kill sign bit
  909. - 
  910. -     andl    #0x0fffff,d4    | remove exponent from u.mantissa
  911. -     andl    #0x0fffff,d6    | remove exponent from v.mantissa
  912.   |
  913.   |
  914.   |
  915. !     cmpw    #0x7ff,d0
  916.       beq    0f        |u == NaN || u== Inf
  917. !     cmpw    #0x7ff,d1
  918.       beq    1f        | v == NaN || v == Inf
  919.       tstw    d0
  920.       bne    3f        | u not zero nor denorm
  921. --- 77,103 ----
  922.       moveml    d2-d7,sp@-    | save registers
  923.       moveml    a0@,d4-d5/d6-d7    | d4-d5 = u, d6-d7 = v
  924.   
  925. +     movel    #0x0fffff,d3
  926.       movel    d4,d0        | d0 = u.exp
  927. +     andl    d3,d4        | remove exponent from u.mantissa
  928.       swap    d0
  929.       movew    d0,d2        | d2 = u.sign
  930. +     bclr    #15,d0        | kill sign bit
  931.       lsrw    #4,d0
  932.   
  933.       movel    d6,d1        | d1 = v.exp
  934. +     andl    d3,d6        | remove exponent from v.mantissa
  935.       swap    d1
  936.       eorw    d1,d2        | d2 = u.sign ^ v.sign (in bit 15)
  937. +     bclr    #15,d1        | kill sign bit
  938.       lsrw    #4,d1
  939.   |
  940.   |
  941.   |
  942. !     movew    #0x7ff,d3
  943. !     cmpw    d3,d0
  944.       beq    0f        |u == NaN || u== Inf
  945. !     cmpw    d3,d1
  946.       beq    1f        | v == NaN || v == Inf
  947.       tstw    d0
  948.       bne    3f        | u not zero nor denorm
  949. ***************
  950. *** 107,117 ****
  951.       orl    d6,d3
  952.       bne    nospec
  953.       bra    retinf        | x/0 -> +/- Inf
  954. !     
  955.   0:    orl    d5,d4        | u == NaN ?
  956.       bne    retnan        | NaN/ x
  957. !     cmpw    #0x7ff,d1    
  958. !     beq    retnan        | Inf/Inf or Inf/NaN 
  959.       bra    retinf        | Inf/x | x != Inf && x != NaN
  960.   
  961.   1:    orl    d7,d6
  962. --- 112,122 ----
  963.       orl    d6,d3
  964.       bne    nospec
  965.       bra    retinf        | x/0 -> +/- Inf
  966. ! 
  967.   0:    orl    d5,d4        | u == NaN ?
  968.       bne    retnan        | NaN/ x
  969. !     cmpw    #0x7ff,d1
  970. !     beq    retnan        | Inf/Inf or Inf/NaN
  971.       bra    retinf        | Inf/x | x != Inf && x != NaN
  972.   
  973.   1:    orl    d7,d6
  974. ***************
  975. *** 121,168 ****
  976.   2:    tstw    d1
  977.       bne    retzero        | 0/x ->+/- 0
  978.       orl    d5,d4
  979. !     bne    retzero        | 0/x 
  980.       bra    retnan        | 0/0
  981.   |
  982.   |    Return Infinity with correct sign
  983. ! |    
  984. ! retinf:    clrl    d1
  985. !     tstw    d2
  986. !     bpl    0f
  987. !     movel    #0xfff00000,d0
  988.   return:    moveml    sp@+,d2-d7
  989.       rts
  990. - 
  991. - 0:    movel    #0x7ff00000,d0
  992. -     bra    return    
  993.   |
  994.   |    Return NaN
  995.   |
  996. ! retnan: movel    #0x7fffffff,d0
  997. !     moveql    #-1,d1
  998.       bra    return
  999.   |
  1000.   |    Return correct signed zero
  1001.   |
  1002. ! retzero:clrl    d0        | zero destination
  1003. !     clrl    d1
  1004. !     tstw    d2
  1005. !     bge    return
  1006. !     bset    #31,d0
  1007.       bra    return
  1008.   |
  1009.   |    End of special handling
  1010. ! |    
  1011. ! nospec:    tstw    d0        | check for zero exponent - no leading "1"
  1012. !     beq    0f
  1013. !     bset    #20,d4        | restore implied leading "1"
  1014. !     bra    1f
  1015. ! 0:    addw    #1,d0        | "normalize" exponent
  1016. ! 
  1017. ! 1:    tstw    d1        | check for zero exponent - no leading "1"
  1018. !     beq    0f
  1019. !     bset    #20,d6        | restore implied leading "1"
  1020. !     bra    1f
  1021.   0:    addw    #1,d1        | "normalize" exponent
  1022.   
  1023.   1:    movew    d2,a0        | save sign
  1024. --- 126,170 ----
  1025.   2:    tstw    d1
  1026.       bne    retzero        | 0/x ->+/- 0
  1027.       orl    d5,d4
  1028. !     bne    retzero        | 0/x
  1029.       bra    retnan        | 0/0
  1030.   |
  1031.   |    Return Infinity with correct sign
  1032. ! |
  1033. ! retinf:    moveq    #0,d1
  1034. !     movel    #0xffe00000,d0
  1035. !     lslw    #1,d2
  1036. !     roxrl   #1,d0        | shift in high bit as given by d2
  1037.   return:    moveml    sp@+,d2-d7
  1038.       rts
  1039.   |
  1040.   |    Return NaN
  1041.   |
  1042. ! retnan: moveql    #-1,d1
  1043. !     movel    d1,d0
  1044. !     bclr    #31,d0        | 0x7fffffff -> d0
  1045.       bra    return
  1046.   |
  1047.   |    Return correct signed zero
  1048.   |
  1049. ! retzero:moveq    #0,d0        | zero destination
  1050. !     movel    d0,d1
  1051. !     lslw    #1,d2        | we need an extension bit
  1052. !     roxrl    #1,d0
  1053.       bra    return
  1054.   |
  1055.   |    End of special handling
  1056. ! |
  1057. ! nospec:    bset    #20,d4        | restore implied leading "1"
  1058. !     tstw    d0        | check for zero exponent - no leading "1"
  1059. !     bne    1f
  1060. !     bclr    #20,d4        | remove it
  1061. !     addw    #1,d0        | "normalize" exponent
  1062. ! 
  1063. ! 1:    bset    #20,d6        | restore implied leading "1"
  1064. !     tstw    d1        | check for zero exponent - no leading "1"
  1065. !     bne    1f
  1066. !     bclr    #20,d6        | remove it
  1067.   0:    addw    #1,d1        | "normalize" exponent
  1068.   
  1069.   1:    movew    d2,a0        | save sign
  1070. ***************
  1071. *** 172,179 ****
  1072.       addw    #66,d0        |  add loop offset, +2 for extra rounding bits
  1073.                   |   for denormalized numbers (2 implied by dbra)
  1074.       movew    #24,d1        | bit number for "implied" pos (+4 for rounding)
  1075. !     movel    #-1,d2        | zero the quotient
  1076. !     movel    #-1,d3        |  (for speed it is a one''s complement)
  1077.       subl    d7,d5        | initial subtraction,
  1078.       subxl    d6,d4        | u = u - v
  1079.   2:
  1080. --- 174,181 ----
  1081.       addw    #66,d0        |  add loop offset, +2 for extra rounding bits
  1082.                   |   for denormalized numbers (2 implied by dbra)
  1083.       movew    #24,d1        | bit number for "implied" pos (+4 for rounding)
  1084. !     moveq    #-1,d2        | zero the quotient
  1085. !     moveq    #-1,d3        |  (for speed it is a one''s complement)
  1086.       subl    d7,d5        | initial subtraction,
  1087.       subxl    d6,d4        | u = u - v
  1088.   2:
  1089. ***************
  1090. *** 196,208 ****
  1091.       addl    d7,d5        | add (restore)
  1092.       addxl    d6,d4        | u = u + v
  1093.       dbra    d0,2b        | give up if result is denormalized
  1094. ! 5:    subw    #2,d0        | remove rounding offset for denormalized nums
  1095.       notl    d2        | invert quotient to get it right
  1096.       notl    d3
  1097.   
  1098.       movel   d5,d1
  1099.       orl     d4,d1           | check for exact result
  1100. !     beq     1f              
  1101.       moveql  #-1,d1          | Set rounding bits for tie case
  1102.   1:    movel    d2,d4        | save quotient mantissa
  1103.       movel    d3,d5
  1104. --- 198,210 ----
  1105.       addl    d7,d5        | add (restore)
  1106.       addxl    d6,d4        | u = u + v
  1107.       dbra    d0,2b        | give up if result is denormalized
  1108. ! 5:    subqw    #2,d0        | remove rounding offset for denormalized nums
  1109.       notl    d2        | invert quotient to get it right
  1110.       notl    d3
  1111.   
  1112.       movel   d5,d1
  1113.       orl     d4,d1           | check for exact result
  1114. !     beq     1f
  1115.       moveql  #-1,d1          | Set rounding bits for tie case
  1116.   1:    movel    d2,d4        | save quotient mantissa
  1117.       movel    d3,d5
  1118. ===================================================================
  1119. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_divmods.s,v
  1120. retrieving revision 1.3
  1121. diff -c -r1.3 _divmods.s
  1122. *** 1.3    1992/11/24 04:14:46
  1123. --- _divmods.s    1993/05/29 20:28:24
  1124. ***************
  1125. *** 62,68 ****
  1126.       cmpl    d2,d0        | compare with divisor
  1127.       bcs    0f
  1128.       subl    d2,d0        | big enough, subtract
  1129. !     addw    #1,d1        | and note bit into result
  1130.   0:
  1131.       dbra    d3,3b
  1132.       exg    d0,d1        | put quotient and remainder in their registers
  1133. --- 62,68 ----
  1134.       cmpl    d2,d0        | compare with divisor
  1135.       bcs    0f
  1136.       subl    d2,d0        | big enough, subtract
  1137. !     addqw    #1,d1        | and note bit into result
  1138.   0:
  1139.       dbra    d3,3b
  1140.       exg    d0,d1        | put quotient and remainder in their registers
  1141. ===================================================================
  1142. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_divsf3.cpp,v
  1143. retrieving revision 1.6
  1144. diff -c -r1.6 _divsf3.cpp
  1145. *** 1.6    1993/05/21 12:55:08
  1146. --- _divsf3.cpp    1993/05/29 20:28:26
  1147. ***************
  1148. *** 18,38 ****
  1149.   ___divsf3:
  1150.   
  1151.   #ifdef    ERROR_CHECK
  1152. !     tstl    a7@(8)            | check if divisor is 0
  1153.       bne    no_exception
  1154.   
  1155.       pea    pc@(LC0)
  1156.       pea    Stderr
  1157. !     jbsr    _fprintf    |
  1158. !     addql    #8,a7        |
  1159. !                     | set _errno to ERANGE
  1160. !     moveq    #ERANGE,d0
  1161.       Emove    d0,Errno
  1162. !     movel    __infinitysf,d0        | return signed infinity 
  1163. !     btst    #31,a7@(4)        | transfer sign of dividend
  1164.       beq    clear            | (mjr++)
  1165.       bset    #31,d0            |
  1166. -     rts                |
  1167.   clear:                    |
  1168.       rts
  1169.   
  1170. --- 18,37 ----
  1171.   ___divsf3:
  1172.   
  1173.   #ifdef    ERROR_CHECK
  1174. !     tstl    sp@(8)            | check if divisor is 0
  1175.       bne    no_exception
  1176.   
  1177.       pea    pc@(LC0)
  1178.       pea    Stderr
  1179. !     jbsr    _fprintf
  1180. !     addql    #8,sp
  1181. ! 
  1182. !     moveq    #Erange,d0        | set _errno to ERANGE
  1183.       Emove    d0,Errno
  1184. !     movel    __infinitysf,d0        | return signed infinity
  1185. !     btst    #31,sp@(4)        | transfer sign of dividend
  1186.       beq    clear            | (mjr++)
  1187.       bset    #31,d0            |
  1188.   clear:                    |
  1189.       rts
  1190.   
  1191. ***************
  1192. *** 43,48 ****
  1193. --- 42,53 ----
  1194.   |
  1195.   | written by Kai-Uwe Bloem (I5110401@dbstu1.bitnet).
  1196.   | Based on a 80x86 floating point packet from comp.os.minix, written by P.Housel
  1197. + | Revision 1.2.4 michal 05-93 (ntomczak@vm.ucs.ualberta.ca)
  1198. + |   + resynchro with errno codes
  1199. + |   + code smoothing
  1200. + |   + removed extra code 'retz:' and 'divz:'; these cases are
  1201. + |     going back through 'retzero:' and 'retinf:'
  1202. + |
  1203.   | patched by Olaf Flebbe (flebbe@tat.physik.uni-tuebingen.de)
  1204.   |
  1205.   | Revision 1.2.3 olaf 5-93
  1206. ***************
  1207. *** 72,103 ****
  1208.       moveml    d2-d5,sp@-    | save registers
  1209.       moveml    a0@,d4/d5    | d4 = u, d5 = v
  1210.   
  1211.       movel    d4,d0        | d0 = u.exp
  1212.       swap    d0
  1213.       movew    d0,d2        | d2 = u.sign
  1214.       lsrw    #7,d0
  1215. -     andw    #0xff,d0    | kill sign bit
  1216.   
  1217.       movel    d5,d1        | d1 = v.exp
  1218.       swap    d1
  1219.       eorw    d1,d2        | d2 = u.sign ^ v.sign (in bit 15)
  1220.       lsrw    #7,d1
  1221. -     andw    #0xff,d1    | kill sign bit
  1222. - 
  1223. -     andl    #0x7fffff,d4    | remove exponent from u.mantissa
  1224. -     andl    #0x7fffff,d5    | remove exponent from v.mantissa
  1225.   |
  1226.   |
  1227.   |
  1228. !     cmpb    #0xff,d0    
  1229.       beq    0f        |u == NaN || u== Inf
  1230. !     cmpb    #0xff,d1
  1231.       beq    1f        | v == NaN || v == Inf
  1232.       tstb    d0
  1233.       bne    3f        | u not zero nor denorm
  1234.       tstl    d4
  1235.       beq    2f        | 0/ ?
  1236. !     
  1237.   3:    tstw    d1
  1238.       bne    nospec
  1239.   
  1240. --- 77,109 ----
  1241.       moveml    d2-d5,sp@-    | save registers
  1242.       moveml    a0@,d4/d5    | d4 = u, d5 = v
  1243.   
  1244. +     movel    #0x7fffff,d3
  1245.       movel    d4,d0        | d0 = u.exp
  1246. +     andl    d3,d4        | remove exponent from u.mantissa
  1247.       swap    d0
  1248.       movew    d0,d2        | d2 = u.sign
  1249. +     bclr    #15,d0        | kill sign bit
  1250.       lsrw    #7,d0
  1251.   
  1252.       movel    d5,d1        | d1 = v.exp
  1253. +     andl    d3,d5        | remove exponent from v.mantissa
  1254.       swap    d1
  1255.       eorw    d1,d2        | d2 = u.sign ^ v.sign (in bit 15)
  1256. +     bclr    #15,d1        | kill sign bit
  1257.       lsrw    #7,d1
  1258.   |
  1259.   |
  1260.   |
  1261. !     moveq    #-1,d3
  1262. !     cmpb    d3,d0        | comparison with #0xff
  1263.       beq    0f        |u == NaN || u== Inf
  1264. !     cmpb    d3,d1
  1265.       beq    1f        | v == NaN || v == Inf
  1266.       tstb    d0
  1267.       bne    3f        | u not zero nor denorm
  1268.       tstl    d4
  1269.       beq    2f        | 0/ ?
  1270. ! 
  1271.   3:    tstw    d1
  1272.       bne    nospec
  1273.   
  1274. ***************
  1275. *** 107,114 ****
  1276.   
  1277.   0:    tstl    d4        | u == NaN ?
  1278.       bne    retnan        | NaN/ x
  1279. !     cmpb    #0xff,d1    
  1280. !     beq    retnan        | Inf/Inf or Inf/NaN 
  1281.       bra    retinf        | Inf/x | x != Inf && x != NaN
  1282.   
  1283.   1:    tstl    d5
  1284. --- 113,120 ----
  1285.   
  1286.   0:    tstl    d4        | u == NaN ?
  1287.       bne    retnan        | NaN/ x
  1288. !     cmpb    d3,d1
  1289. !     beq    retnan        | Inf/Inf or Inf/NaN
  1290.       bra    retinf        | Inf/x | x != Inf && x != NaN
  1291.   
  1292.   1:    tstl    d5
  1293. ***************
  1294. *** 118,177 ****
  1295.   2:    tstw    d1
  1296.       bne    retzero        | 0/x ->+/- 0
  1297.       tstl    d4
  1298. !     bne    retzero        | 0/x 
  1299.       bra    retnan        | 0/0
  1300.   |
  1301.   |    Return Infinity with correct sign
  1302. ! |    
  1303. ! retinf:    tstw    d2
  1304. !     bpl    0f
  1305. !     movel    #0xff800000,d0
  1306.   return:    moveml    sp@+,d2-d5
  1307.       rts
  1308. - 
  1309. - 0:    movel    #0x7f800000,d0
  1310. -     bra    return    
  1311.   |
  1312.   |    Return NaN
  1313.   |
  1314. ! retnan: movel    #0x7fffffff,d0
  1315.       bra    return
  1316.   |
  1317.   |    Return correct signed zero
  1318.   |
  1319. ! retzero:clrl    d0        | zero destination
  1320. !     tstw    d2
  1321. !     bge    return
  1322. !     bset    #31,d0
  1323.       bra    return
  1324.   |
  1325.   |    End of special handling
  1326. ! |    
  1327. ! nospec:    tstw    d0        | check for zero exponent - no leading "1"
  1328. !     beq    0f
  1329. !     bset    #23,d4        | restore implied leading "1"
  1330. !     bra    1f
  1331. ! 0:    addw    #1,d0        | "normalize" exponent
  1332.   1:
  1333. - # ifndef ERROR_CHECK
  1334.       tstl    d4
  1335. !     beq    retz        | dividing zero
  1336. ! # endif    ERROR_CHECK
  1337.   
  1338. -     tstw    d1        | check for zero exponent - no leading "1"
  1339. -     beq    0f
  1340.       bset    #23,d5        | restore implied leading "1"
  1341. !     bra    1f
  1342. ! 0:    addw    #1,d1        | "normalize" exponent
  1343. ! 1:    tstl    d5
  1344. !     beq    divz        | divide by zero
  1345.   
  1346.       subw    d1,d0        | subtract exponents,
  1347.       addw    #BIAS4-8+1,d0    |  add bias back in, account for shift
  1348.       addw    #34,d0        |  add loop offset, +2 for extra rounding bits
  1349.                   |   for denormalized numbers (2 implied by dbra)
  1350.       movew    #27,d1        | bit number for "implied" pos (+4 for rounding)
  1351. !     movel    #-1,d3        |  zero quotient (for speed a one''s complement)
  1352.       subl    d5,d4        | initial subtraction, u = u - v
  1353.   2:
  1354.       btst    d1,d3        | divide until 1 in implied position
  1355. --- 124,181 ----
  1356.   2:    tstw    d1
  1357.       bne    retzero        | 0/x ->+/- 0
  1358.       tstl    d4
  1359. !     bne    retzero        | 0/x
  1360.       bra    retnan        | 0/0
  1361.   |
  1362.   |    Return Infinity with correct sign
  1363. ! |
  1364. ! retinf:    movel    #0xff000000,d0
  1365. !     lslw    #1,d2
  1366. !     roxrl   #1,d0        | shift in high bit as given by d2
  1367.   return:    moveml    sp@+,d2-d5
  1368.       rts
  1369.   |
  1370.   |    Return NaN
  1371.   |
  1372. ! retnan: movel    d3,d0        | d3 contains 0xffffffff
  1373. !     bclr    #31,d0
  1374.       bra    return
  1375.   |
  1376.   |    Return correct signed zero
  1377.   |
  1378. ! retzero:moveq    #0,d0        | zero destination
  1379. !     lslw    #1,d2        | set X bit accordingly
  1380. !     roxrl    #1,d0
  1381.       bra    return
  1382.   |
  1383.   |    End of special handling
  1384. ! |
  1385. ! nospec:    bset    #23,d4        | restore implied leading "1"
  1386. !     tstw    d0        | check for zero exponent - no leading "1"
  1387. !     bne    1f
  1388. !     bclr    #23,d4        | remove it
  1389. !     addw    #1,d0        | "normalize" exponent
  1390.   1:
  1391.       tstl    d4
  1392. !     beq    retzero        | dividing zero
  1393.   
  1394.       bset    #23,d5        | restore implied leading "1"
  1395. !     tstw    d1        | check for zero exponent - no leading "1"
  1396. !     bne    1f
  1397. !     bclr    #23,d5        | remove it
  1398. !     addw    #1,d1        | "normalize" exponent
  1399. ! 1:
  1400. ! # ifndef ERROR_CHECK
  1401. !     tstl    d5
  1402. !     beq    retinf        | divide by zero
  1403. ! # endif    ERROR_CHECK
  1404.   
  1405.       subw    d1,d0        | subtract exponents,
  1406.       addw    #BIAS4-8+1,d0    |  add bias back in, account for shift
  1407.       addw    #34,d0        |  add loop offset, +2 for extra rounding bits
  1408.                   |   for denormalized numbers (2 implied by dbra)
  1409.       movew    #27,d1        | bit number for "implied" pos (+4 for rounding)
  1410. !     moveql    #-1,d3        |  zero quotient (for speed a one''s complement)
  1411.       subl    d5,d4        | initial subtraction, u = u - v
  1412.   2:
  1413.       btst    d1,d3        | divide until 1 in implied position
  1414. ***************
  1415. *** 188,194 ****
  1416.       addxl    d3,d3        | shift quotient and clear bit zero
  1417.       addl    d5,d4        | add (restore), u = u + v
  1418.       dbra    d0,2b        | give up if result is denormalized
  1419. ! 5:    subw    #2,d0        | remove rounding offset for denormalized nums
  1420.       notl    d3        | invert quotient to get it right
  1421.   
  1422.       clrl    d1        | zero rounding bits
  1423. --- 192,198 ----
  1424.       addxl    d3,d3        | shift quotient and clear bit zero
  1425.       addl    d5,d4        | add (restore), u = u + v
  1426.       dbra    d0,2b        | give up if result is denormalized
  1427. ! 5:    subqw    #2,d0        | remove rounding offset for denormalized nums
  1428.       notl    d3        | invert quotient to get it right
  1429.   
  1430.       clrl    d1        | zero rounding bits
  1431. ***************
  1432. *** 198,205 ****
  1433.   1:    movel    d3,d4        | save quotient mantissa
  1434.       jmp    norm_sf        | (registers on stack removed by norm_sf)
  1435.   
  1436.   # ifndef ERROR_CHECK
  1437. ! retz:    clrl    d0        | zero destination
  1438.       moveml    sp@+,d2-d5
  1439.       rts            | no normalization needed
  1440.   
  1441. --- 202,211 ----
  1442.   1:    movel    d3,d4        | save quotient mantissa
  1443.       jmp    norm_sf        | (registers on stack removed by norm_sf)
  1444.   
  1445. + #if 0
  1446. + | this  is dead code right now - retzero and retinf used (mj)
  1447.   # ifndef ERROR_CHECK
  1448. ! retz:    moveq    #0,d0        | zero destination
  1449.       moveml    sp@+,d2-d5
  1450.       rts            | no normalization needed
  1451.   
  1452. ***************
  1453. *** 213,222 ****
  1454.       bclr    #31,d0        |
  1455.       rts
  1456.   # endif    ERROR_CHECK
  1457.   #else
  1458.   
  1459.   | single precision floating point stuff for Atari-gcc using the SFP004
  1460. ! | or compatible boards with a memory mapped 68881 
  1461.   | developed with gas
  1462.   |
  1463.   |  single floating point divide routine
  1464. --- 219,230 ----
  1465.       bclr    #31,d0        |
  1466.       rts
  1467.   # endif    ERROR_CHECK
  1468. + #endif
  1469. + 
  1470.   #else
  1471.   
  1472.   | single precision floating point stuff for Atari-gcc using the SFP004
  1473. ! | or compatible boards with a memory mapped 68881
  1474.   | developed with gas
  1475.   |
  1476.   |  single floating point divide routine
  1477. ***************
  1478. *** 247,256 ****
  1479.       lea    0xfffffa50:w,a0
  1480.       movew    #0x4400,a0@(comm)    | load first argument to fp0
  1481.       cmpiw    #0x8900,a0@(resp)    | check
  1482. !     movel    a7@(4),a0@
  1483.       movew    #0x4424,a0@(comm)
  1484.       .long    0x0c688900, 0xfff067f8
  1485. !     movel    a7@(8),a0@
  1486.       movew    #0x6400,a0@(comm)    | result to d0
  1487.       .long    0x0c688900, 0xfff067f8
  1488.       movel    a0@,d0            | REMARK: 0/0 returns a NAN
  1489. --- 255,264 ----
  1490.       lea    0xfffffa50:w,a0
  1491.       movew    #0x4400,a0@(comm)    | load first argument to fp0
  1492.       cmpiw    #0x8900,a0@(resp)    | check
  1493. !     movel    sp@(4),a0@
  1494.       movew    #0x4424,a0@(comm)
  1495.       .long    0x0c688900, 0xfff067f8
  1496. !     movel    sp@(8),a0@
  1497.       movew    #0x6400,a0@(comm)    | result to d0
  1498.       .long    0x0c688900, 0xfff067f8
  1499.       movel    a0@,d0            | REMARK: 0/0 returns a NAN
  1500. ===================================================================
  1501. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_extends.cpp,v
  1502. retrieving revision 1.4
  1503. diff -c -r1.4 _extends.cpp
  1504. *** 1.4    1993/03/01 17:34:03
  1505. --- _extends.cpp    1993/05/29 20:28:27
  1506. ***************
  1507. *** 14,19 ****
  1508. --- 14,22 ----
  1509.   | written by Kai-Uwe Bloem (I5110401@dbstu1.bitnet).
  1510.   | Based on a 80x86 floating point packet from comp.os.minix, written by P.Housel
  1511.   |
  1512. + | Revision 1.3, michal 05-93 :
  1513. + | code cleanup
  1514. + |
  1515.   | Revision 1.2, kub 01-90 :
  1516.   | added support for denormalized numbers
  1517.   |
  1518. ***************
  1519. *** 26,64 ****
  1520.   BIAS4    =    0x7F-1
  1521.   BIAS8    =    0x3FF-1
  1522.   
  1523. !     lea    sp@(4),a0    | parameter pointer
  1524.       moveml    d2-d7,sp@-    | save regs to keep norm_df happy
  1525. !     movel    a0@,d4        | get number
  1526. !     clrl    d5        | prepare double mantissa
  1527.   
  1528. !     movew    a0@,d0        | extract exponent
  1529.       movew    d0,d2        | extract sign
  1530.       lsrw    #7,d0
  1531. -     andw    #0xff,d0    | kill sign bit (exponent is 8 bits)
  1532.       andl    #0x7fffff,d4    | remove exponent from mantissa
  1533.   
  1534. !     cmpb    #0xff,d0
  1535.       bne    nospec
  1536.       tstl    d4
  1537.       beq    retinf
  1538.   | ret nan
  1539. !     movel    #0x7fffffff,d0
  1540. !     moveql    #-1,d1
  1541.   return:    moveml    sp@+,d2-d7
  1542.       rts
  1543. - retinf:    clrl    d1
  1544. -     movel    #0x7ff00000,d0
  1545. -     tstw    d2
  1546. -     bpl    return
  1547. -     bset    #31,d0
  1548. -     bra    return
  1549.   
  1550. !     
  1551. ! nospec:    tstw    d0        | check for zero exponent - no leading "1"
  1552. !     beq    0f        | for denormalized numbers
  1553. !     bset    #23,d4        | restore implied leading "1"
  1554. !     bra    1f
  1555. ! 0:    addw    #1,d0        | "normalize" exponent
  1556.   1:
  1557.       addw    #BIAS8-BIAS4-3,d0    | adjust bias, account for shift
  1558.       clrw    d1        | dummy rounding info
  1559. --- 29,66 ----
  1560.   BIAS4    =    0x7F-1
  1561.   BIAS8    =    0x3FF-1
  1562.   
  1563. !     movel   sp@(4),d0    | get number
  1564.       moveml    d2-d7,sp@-    | save regs to keep norm_df happy
  1565. !     moveq    #0,d5        | prepare double mantissa
  1566.   
  1567. !     movel    d0,d4
  1568. !     swap    d0        | extract exponent
  1569.       movew    d0,d2        | extract sign
  1570. +     bclr    #15,d0        | kill sign bit (exponent is 8 bits)
  1571.       lsrw    #7,d0
  1572.       andl    #0x7fffff,d4    | remove exponent from mantissa
  1573.   
  1574. !     moveql    #-1,d1
  1575. !     cmpb    d1,d0
  1576.       bne    nospec
  1577.       tstl    d4
  1578.       beq    retinf
  1579.   | ret nan
  1580. !     movel    d1,d0
  1581. !     lsrl    #1,d0        | #0x7fffffff -> d0
  1582. !     bra    return
  1583. ! retinf:    movel    d4,d1        | 0 -> d1
  1584. !     movel    #0xffe00000,d0
  1585. !     lslw    #1,d2        | get extension bit
  1586. !     roxrl    #1,d0        | shift in sign bit
  1587.   return:    moveml    sp@+,d2-d7
  1588.       rts
  1589.   
  1590. ! nospec: bset    #23,d4        | restore implied leading "1"
  1591. !     tstw    d0        | check for zero exponent - no leading "1"
  1592. !     bne    1f        | for denormalized numbers
  1593. !     bclr    #23,d4        | so clear it and ...
  1594. !     addw    #1,d0        | "normalize" exponent
  1595.   1:
  1596.       addw    #BIAS8-BIAS4-3,d0    | adjust bias, account for shift
  1597.       clrw    d1        | dummy rounding info
  1598. ***************
  1599. *** 97,103 ****
  1600.       lea    0xfffffa50:w,a0
  1601.       movew    #0x4400,a0@(comm)    | load argument to fp0
  1602.       cmpiw    #0x8900,a0@(resp)    | check
  1603. !     movel    a7@(4),a0@        | now push arg
  1604.       movew    #0x7400,a0@(comm)    | result to d0/d1
  1605.       .long    0x0c688900, 0xfff067f8
  1606.       movel    a0@,d0            | pop double float
  1607. --- 99,105 ----
  1608.       lea    0xfffffa50:w,a0
  1609.       movew    #0x4400,a0@(comm)    | load argument to fp0
  1610.       cmpiw    #0x8900,a0@(resp)    | check
  1611. !     movel    sp@(4),a0@        | now push arg
  1612.       movew    #0x7400,a0@(comm)    | result to d0/d1
  1613.       .long    0x0c688900, 0xfff067f8
  1614.       movel    a0@,d0            | pop double float
  1615. ***************
  1616. *** 122,130 ****
  1617.   | no NAN checking implemented since the 68881 treats this situation "correct",
  1618.   | i.e. according to IEEE
  1619.   
  1620. !     fmoves    a7@(4),fp0    | load argument to fp0
  1621. !     fmoved    fp0,a7@-    | read back as double
  1622. !     moveml    a7@+,d0-d1
  1623.       rts
  1624.   
  1625.   #endif    __M68881__
  1626. --- 124,132 ----
  1627.   | no NAN checking implemented since the 68881 treats this situation "correct",
  1628.   | i.e. according to IEEE
  1629.   
  1630. !     fmoves    sp@(4),fp0    | load argument to fp0
  1631. !     fmoved    fp0,sp@-    | read back as double
  1632. !     moveml    sp@+,d0-d1
  1633.       rts
  1634.   
  1635.   #endif    __M68881__
  1636. ===================================================================
  1637. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_fixdfsi.cpp,v
  1638. retrieving revision 1.4
  1639. diff -c -r1.4 _fixdfsi.cpp
  1640. *** 1.4    1993/03/01 17:34:05
  1641. --- _fixdfsi.cpp    1993/05/29 20:28:28
  1642. ***************
  1643. *** 40,51 ****
  1644.       lea    0xfffffa50:w,a0
  1645.       movew    #0x5403,a0@(comm)    | fintrz to fp0
  1646.       cmpiw    #0x8900,a0@(resp)    | check
  1647. !     movel    a7@(4),a0@
  1648. !     movel    a7@(8),a0@
  1649.       movew    #0x6000,a0@(comm)    | result to d0
  1650.       .long    0x0c688900, 0xfff067f8
  1651.       movel    a0@,d0
  1652. !         
  1653.   #endif    sfp004
  1654.   #ifdef    __M68881__
  1655.   
  1656. --- 40,51 ----
  1657.       lea    0xfffffa50:w,a0
  1658.       movew    #0x5403,a0@(comm)    | fintrz to fp0
  1659.       cmpiw    #0x8900,a0@(resp)    | check
  1660. !     movel    sp@(4),a0@
  1661. !     movel    sp@(8),a0@
  1662.       movew    #0x6000,a0@(comm)    | result to d0
  1663.       .long    0x0c688900, 0xfff067f8
  1664.       movel    a0@,d0
  1665. ! 
  1666.   #endif    sfp004
  1667.   #ifdef    __M68881__
  1668.   
  1669. ***************
  1670. *** 98,111 ****
  1671.       moveml    a0@,d0-d1    | get the number
  1672.       movew    a0@,d2        | extract exp
  1673.       movew    d2,d3        | extract sign
  1674.       lsrw    #4,d2
  1675. -     andw    #0x07ff,d2    | kill sign bit
  1676.   
  1677.       andl    #0x0fffff,d0    | remove exponent from mantissa
  1678.       bset    #20,d0        | restore implied leading "1"
  1679.   
  1680.       cmpw    #BIAS8,d2    | check exponent
  1681. !     blt    zero        | strictly factional, no integer part ?
  1682.       cmpw    #BIAS8+32,d2    | is it too big to fit in a 32-bit integer ?
  1683.       bgt    toobig
  1684.   
  1685. --- 98,111 ----
  1686.       moveml    a0@,d0-d1    | get the number
  1687.       movew    a0@,d2        | extract exp
  1688.       movew    d2,d3        | extract sign
  1689. +     bclr    #15,d2        | kill sign bit
  1690.       lsrw    #4,d2
  1691.   
  1692.       andl    #0x0fffff,d0    | remove exponent from mantissa
  1693.       bset    #20,d0        | restore implied leading "1"
  1694.   
  1695.       cmpw    #BIAS8,d2    | check exponent
  1696. !     blt    zero        | strictly fractional, no integer part ?
  1697.       cmpw    #BIAS8+32,d2    | is it too big to fit in a 32-bit integer ?
  1698.       bgt    toobig
  1699.   
  1700. ***************
  1701. *** 120,126 ****
  1702.   
  1703.   2:    addl    d1,d1        | shift up to align radix point
  1704.       addxl    d0,d0
  1705. !     subw    #1,d2
  1706.       bgt    2b
  1707.   
  1708.   3:    cmpl    #0x80000000,d0    | -2147483648 is a nasty evil special case
  1709. --- 120,126 ----
  1710.   
  1711.   2:    addl    d1,d1        | shift up to align radix point
  1712.       addxl    d0,d0
  1713. !     subqw    #1,d2
  1714.       bgt    2b
  1715.   
  1716.   3:    cmpl    #0x80000000,d0    | -2147483648 is a nasty evil special case
  1717. ***************
  1718. *** 143,149 ****
  1719.       bra    8b
  1720.   
  1721.   toobig:
  1722. !     movel    #0x7fffffff,d0    | ugh. Should cause a trap here.
  1723.       moveml    sp@+,d2/d3
  1724.   
  1725.   #endif /* !defined (sfp004) && !defined (__M68881__) */
  1726. --- 143,150 ----
  1727.       bra    8b
  1728.   
  1729.   toobig:
  1730. !     moveq    #-1,d0        | ugh. Should cause a trap here.
  1731. !     bclr    #31,d0
  1732.       moveml    sp@+,d2/d3
  1733.   
  1734.   #endif /* !defined (sfp004) && !defined (__M68881__) */
  1735. ***************
  1736. *** 157,171 ****
  1737.       ble    error_msg
  1738.       rts
  1739.   error_msg:
  1740. !     moveml    d0-d1,a7@-
  1741. !     moveq    #ERANGE,d0
  1742.       Emove    d0,Errno
  1743.       pea    pc@(_Overflow)    | for printf
  1744. !     pea    pc@(_Error_String)    |
  1745.       pea    Stderr
  1746. !     jbsr    _fprintf    |
  1747. !     addl    #12,a7        |
  1748. !     moveml    a7@+,d0-d1
  1749.   #endif    /* ERROR_CHECK */
  1750.   
  1751.       rts
  1752. --- 158,172 ----
  1753.       ble    error_msg
  1754.       rts
  1755.   error_msg:
  1756. !     moveml    d0-d1,sp@-
  1757. !     moveq    #Erange,d0
  1758.       Emove    d0,Errno
  1759.       pea    pc@(_Overflow)    | for printf
  1760. !     pea    pc@(_Error_String)
  1761.       pea    Stderr
  1762. !     jbsr    _fprintf
  1763. !     addl    #12,sp
  1764. !     moveml    sp@+,d0-d1
  1765.   #endif    /* ERROR_CHECK */
  1766.   
  1767.       rts
  1768. ===================================================================
  1769. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_fixsfsi.cpp,v
  1770. retrieving revision 1.3
  1771. diff -c -r1.3 _fixsfsi.cpp
  1772. *** 1.3    1993/03/01 17:34:07
  1773. --- _fixsfsi.cpp    1993/05/29 20:28:30
  1774. ***************
  1775. *** 28,34 ****
  1776.   1:    cmpiw    #0x8900,a0@(resp)
  1777.       beq    1b
  1778.       movel    a0@,d0
  1779. !         
  1780.   #else /* !sfp004 */
  1781.   
  1782.   BIAS4    =    0x7F-1
  1783. --- 28,34 ----
  1784.   1:    cmpiw    #0x8900,a0@(resp)
  1785.       beq    1b
  1786.       movel    a0@,d0
  1787. ! 
  1788.   #else /* !sfp004 */
  1789.   
  1790.   BIAS4    =    0x7F-1
  1791. ***************
  1792. *** 38,45 ****
  1793.       movel    d0,d1
  1794.       swap    d1        | extract exp
  1795.       movew    d1,d2        | extract sign
  1796.       lsrw    #7,d1
  1797. -     andw    #0xff,d1    | kill sign bit
  1798.   
  1799.       andl    #0x7fffff,d0    | remove exponent from mantissa
  1800.       bset    #23,d0        | restore implied leading "1"
  1801. --- 38,45 ----
  1802.       movel    d0,d1
  1803.       swap    d1        | extract exp
  1804.       movew    d1,d2        | extract sign
  1805. +     bclr    #15,d1        | kill sign bit
  1806.       lsrw    #7,d1
  1807.   
  1808.       andl    #0x7fffff,d0    | remove exponent from mantissa
  1809.       bset    #23,d0        | restore implied leading "1"
  1810. ***************
  1811. *** 80,86 ****
  1812.       bra    8b
  1813.   
  1814.   toobig:
  1815. !     movel    #0x7fffffff,d0    | ugh. Should cause a trap here.
  1816.       bra    8b
  1817.   
  1818.   #endif /* !sfp004*/
  1819. --- 80,87 ----
  1820.       bra    8b
  1821.   
  1822.   toobig:
  1823. !     moveq    #-1,d0        | ugh. Should cause a trap here.
  1824. !     bclr    #31,d0        | make it #0x7fffffff
  1825.       bra    8b
  1826.   
  1827.   #endif /* !sfp004*/
  1828. ===================================================================
  1829. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_fltsi.cpp,v
  1830. retrieving revision 1.3
  1831. diff -c -r1.3 _fltsi.cpp
  1832. *** 1.3    1992/10/09 20:35:29
  1833. --- _fltsi.cpp    1993/05/29 20:28:31
  1834. ***************
  1835. *** 72,86 ****
  1836.       moveml    d2-d7,sp@-    | save registers to make norm_df happy
  1837.   
  1838.       movel    d0,d4        | prepare result mantissa
  1839. !     clrl    d5
  1840.       movew    #BIAS8+32-11,d0    | radix point after 32 bits
  1841.   0:
  1842.       movel    d4,d2        | set sign flag
  1843. -     swap    d2
  1844. -     tstw    d2        | check sign of number
  1845.       bge    1f        | nonnegative
  1846.       negl    d4        | take absolute value
  1847.   1:
  1848.       clrw    d1        | set rounding = 0
  1849.       jmp    norm_df
  1850.   
  1851. --- 72,85 ----
  1852.       moveml    d2-d7,sp@-    | save registers to make norm_df happy
  1853.   
  1854.       movel    d0,d4        | prepare result mantissa
  1855. !     moveq    #0,d5
  1856.       movew    #BIAS8+32-11,d0    | radix point after 32 bits
  1857.   0:
  1858.       movel    d4,d2        | set sign flag
  1859.       bge    1f        | nonnegative
  1860.       negl    d4        | take absolute value
  1861.   1:
  1862. +     swap    d2        | follow norm_df conventions
  1863.       clrw    d1        | set rounding = 0
  1864.       jmp    norm_df
  1865.   
  1866. ===================================================================
  1867. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_fltsisf.cpp,v
  1868. retrieving revision 1.2
  1869. diff -c -r1.2 _fltsisf.cpp
  1870. *** 1.2    1992/10/09 20:35:29
  1871. --- _fltsisf.cpp    1993/05/29 20:28:32
  1872. ***************
  1873. *** 39,49 ****
  1874.       movel    sp@(20),d4    | prepare result mantissa
  1875.       movew    #BIAS4+32-8,d0    | radix point after 32 bits
  1876.       movel    d4,d2        | set sign flag
  1877. -     swap    d2
  1878. -     tstw    d2        | check sign of number
  1879.       bge    1f        | nonnegative
  1880.       negl    d4        | take absolute value
  1881.   1:
  1882.       clrw    d1        | set rounding = 0
  1883.       jmp    norm_sf
  1884.   
  1885. --- 39,48 ----
  1886.       movel    sp@(20),d4    | prepare result mantissa
  1887.       movew    #BIAS4+32-8,d0    | radix point after 32 bits
  1888.       movel    d4,d2        | set sign flag
  1889.       bge    1f        | nonnegative
  1890.       negl    d4        | take absolute value
  1891.   1:
  1892. +     swap    d2        | follow norm_sf conventions
  1893.       clrw    d1        | set rounding = 0
  1894.       jmp    norm_sf
  1895.   
  1896. ===================================================================
  1897. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_fxunsd.cpp,v
  1898. retrieving revision 1.3
  1899. diff -c -r1.3 _fxunsd.cpp
  1900. *** 1.3    1993/03/01 17:34:09
  1901. --- _fxunsd.cpp    1993/05/29 20:28:33
  1902. ***************
  1903. *** 21,29 ****
  1904.   ___fixunsdfsi:
  1905.   
  1906.   #ifdef ERROR_CHECK
  1907. !     tstl    a7@(4)            | negative?
  1908.       beq    Continue
  1909. !     moveq    #ERANGE,d0
  1910.       Emove    d0,Errno
  1911.       pea    pc@(_Negative)
  1912.       bra    error_exit
  1913. --- 21,29 ----
  1914.   ___fixunsdfsi:
  1915.   
  1916.   #ifdef ERROR_CHECK
  1917. !     tstl    sp@(4)            | negative?
  1918.       beq    Continue
  1919. !     moveq    #Erange,d0
  1920.       Emove    d0,Errno
  1921.       pea    pc@(_Negative)
  1922.       bra    error_exit
  1923. ***************
  1924. *** 31,37 ****
  1925.   #endif /* ERROR_CHECK */
  1926.   #ifdef __M68881__
  1927.   
  1928. !     fintrzd a7@(4),fp0        | convert
  1929.       fmovel    fp0,d0
  1930.   
  1931.   #endif __M68881__
  1932. --- 31,37 ----
  1933.   #endif /* ERROR_CHECK */
  1934.   #ifdef __M68881__
  1935.   
  1936. !     fintrzd sp@(4),fp0        | convert
  1937.       fmovel    fp0,d0
  1938.   
  1939.   #endif __M68881__
  1940. ***************
  1941. *** 48,55 ****
  1942.       lea    0xfffffa50:w,a0
  1943.       movew    #0x5403,a0@(comm)    | fintrz to fp0
  1944.       cmpiw    #0x8900,a0@(resp)    | check
  1945. !     movel    a7@(4),a0@
  1946. !     movel    a7@(8),a0@
  1947.       movew    #0x6000,a0@(comm)    | result to d0
  1948.       .long    0x0c688900, 0xfff067f8
  1949.       movel    a0@,d0
  1950. --- 48,55 ----
  1951.       lea    0xfffffa50:w,a0
  1952.       movew    #0x5403,a0@(comm)    | fintrz to fp0
  1953.       cmpiw    #0x8900,a0@(resp)    | check
  1954. !     movel    sp@(4),a0@
  1955. !     movel    sp@(8),a0@
  1956.       movew    #0x6000,a0@(comm)    | result to d0
  1957.       .long    0x0c688900, 0xfff067f8
  1958.       movel    a0@,d0
  1959. ***************
  1960. *** 81,94 ****
  1961.       moveml    d4/d5,sp@-    | save registers
  1962.       moveml    a0@,d4-d5    | get the number
  1963.       movew    a0@,d0        | extract exp
  1964.       lsrw    #4,d0
  1965. -     andw    #0x07ff,d0    | kill sign bit
  1966.   
  1967.       andl    #0x0fffff,d4    | remove exponent from mantissa
  1968.       bset    #20,d4        | restore implied leading "1"
  1969.   
  1970.       cmpw    #BIAS8,d0    | check exponent
  1971. !     blt    zero        | strictly factional, no integer part ?
  1972.       cmpw    #BIAS8+32,d0    | is it too big to fit in a 32-bit integer ?
  1973.       bgt    toobig
  1974.   
  1975. --- 81,94 ----
  1976.       moveml    d4/d5,sp@-    | save registers
  1977.       moveml    a0@,d4-d5    | get the number
  1978.       movew    a0@,d0        | extract exp
  1979. +     bclr    #15,d0        | kill sign bit
  1980.       lsrw    #4,d0
  1981.   
  1982.       andl    #0x0fffff,d4    | remove exponent from mantissa
  1983.       bset    #20,d4        | restore implied leading "1"
  1984.   
  1985.       cmpw    #BIAS8,d0    | check exponent
  1986. !     blt    zero        | strictly fractional, no integer part ?
  1987.       cmpw    #BIAS8+32,d0    | is it too big to fit in a 32-bit integer ?
  1988.       bgt    toobig
  1989.   
  1990. ***************
  1991. *** 107,119 ****
  1992.       beq    3f
  1993.   
  1994.   1:    lsrl    #1,d4        | shift down to align radix point;
  1995. !     addw    #1,d0        | extra bits fall off the end (no rounding)
  1996.       blt    1b        | shifted all the way down yet ?
  1997.       bra    3f
  1998.   
  1999.   2:    addl    d5,d5        | shift up to align radix point
  2000.       addxl    d4,d4
  2001. !     subw    #1,d0
  2002.       bgt    2b
  2003.   3:
  2004.       movel    d4,d0        | put integer into result register
  2005. --- 107,119 ----
  2006.       beq    3f
  2007.   
  2008.   1:    lsrl    #1,d4        | shift down to align radix point;
  2009. !     addqw    #1,d0        | extra bits fall off the end (no rounding)
  2010.       blt    1b        | shifted all the way down yet ?
  2011.       bra    3f
  2012.   
  2013.   2:    addl    d5,d5        | shift up to align radix point
  2014.       addxl    d4,d4
  2015. !     subqw    #1,d0
  2016.       bgt    2b
  2017.   3:
  2018.       movel    d4,d0        | put integer into result register
  2019. ***************
  2020. *** 127,133 ****
  2021.   
  2022.   toobig:
  2023.       moveml    sp@+,d4/d5
  2024. !     movel    #0x7fffffff,d0    | ugh. Should cause a trap here.
  2025.   #endif
  2026.   
  2027.   #ifdef    ERROR_CHECK
  2028. --- 127,134 ----
  2029.   
  2030.   toobig:
  2031.       moveml    sp@+,d4/d5
  2032. !     moveq    #-1,d0        | ugh. Should cause a trap here.
  2033. !     bclr    #31,d0
  2034.   #endif
  2035.   
  2036.   #ifdef    ERROR_CHECK
  2037. ***************
  2038. *** 135,150 ****
  2039.       bge    error_plus    |
  2040.       rts
  2041.   error_plus:
  2042. !     moveml    d0-d1,a7@-
  2043. !     moveq    #ERANGE,d0
  2044.       Emove    d0,Errno
  2045.       pea    pc@(_Overflow)    | for printf
  2046.   error_exit:
  2047.       pea    pc@(_Error_String)    |
  2048.       pea    Stderr
  2049.       jbsr    _fprintf    |
  2050. !     addl    #12,a7        |
  2051. !     moveml    a7@+,d0-d1
  2052.   #endif    ERROR_CHECK
  2053.       rts
  2054.   
  2055. --- 136,151 ----
  2056.       bge    error_plus    |
  2057.       rts
  2058.   error_plus:
  2059. !     moveml    d0-d1,sp@-
  2060. !     moveq    #Erange,d0
  2061.       Emove    d0,Errno
  2062.       pea    pc@(_Overflow)    | for printf
  2063.   error_exit:
  2064.       pea    pc@(_Error_String)    |
  2065.       pea    Stderr
  2066.       jbsr    _fprintf    |
  2067. !     addl    #12,sp        |
  2068. !     moveml    sp@+,d0-d1
  2069.   #endif    ERROR_CHECK
  2070.       rts
  2071.   
  2072. ===================================================================
  2073. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_muldf3.cpp,v
  2074. retrieving revision 1.4
  2075. diff -c -r1.4 _muldf3.cpp
  2076. *** 1.4    1993/03/01 17:34:12
  2077. --- _muldf3.cpp    1993/05/29 20:28:35
  2078. ***************
  2079. *** 58,65 ****
  2080. --- 58,72 ----
  2081.   |
  2082.   | written by Kai-Uwe Bloem (I5110401@dbstu1.bitnet).
  2083.   | Based on a 80x86 floating point packet from comp.os.minix, written by P.Housel
  2084. + |
  2085. + | Revision 1.2.3 michal 05-93 (ntomczak@vm.ucs.ualberta.ca)
  2086. + |   + code smoothing
  2087. + |
  2088.   | patched by Olaf Flebbe (flebbe@tat.physik.uni-tuebingen.de)
  2089.   |
  2090. + | Revision 1.2.2 olaf 05-93:
  2091. + |   + fixed a bug for signed bug for 0.
  2092. + |
  2093.   | Revision 1.2.1 olaf 12-92:
  2094.   |   + added support for NaN and Infinites
  2095.   |   + added support for -0
  2096. ***************
  2097. *** 79,104 ****
  2098.       moveml    d2-d7,sp@-
  2099.       moveml    a0@,d4-d5/d6-d7 | d4-d5 = v, d6-d7 = u
  2100.   
  2101.       movel    d6,d0        | d0 = u.exp
  2102.       swap    d0
  2103.       movew    d0,d2        | d2 = u.sign
  2104.       lsrw    #4,d0
  2105. -     andw    #0x07ff,d0    | kill sign bit
  2106.   
  2107.       movel    d4,d1        | d1 = v.exp
  2108.       swap    d1
  2109. !     eorw    d1,d2        | d2 = u.sign ^ v.sign (in bit 31)
  2110.       lsrw    #4,d1
  2111. -     andw    #0x07ff,d1    | kill sign bit
  2112. - 
  2113. -     andl    #0x0fffff,d6    | remove exponent from u.mantissa
  2114. -     andl    #0x0fffff,d4    | remove exponent from v.mantissa
  2115.   |
  2116.   | Testing for NaN and Infinities
  2117.   |
  2118. !     cmpw    #0x7ff,d0
  2119.       beq    0f
  2120. !     cmpw    #0x7ff,d1
  2121.       bne    nospec
  2122.       bra    1f
  2123.   |    first operand is special
  2124. --- 86,112 ----
  2125.       moveml    d2-d7,sp@-
  2126.       moveml    a0@,d4-d5/d6-d7 | d4-d5 = v, d6-d7 = u
  2127.   
  2128. +     movel    #0x0fffff,d3
  2129.       movel    d6,d0        | d0 = u.exp
  2130. +     andl    d3,d6        | remove exponent from u.mantissa
  2131.       swap    d0
  2132.       movew    d0,d2        | d2 = u.sign
  2133. +     bclr    #15,d0        | kill sign bit
  2134.       lsrw    #4,d0
  2135.   
  2136.       movel    d4,d1        | d1 = v.exp
  2137. +     andl    d3,d4        | remove exponent from v.mantissa
  2138.       swap    d1
  2139. !     eorw    d1,d2        | d2 = u.sign ^ v.sign (in bit 15)
  2140. !     bclr    #15,d1        | kill sign bit
  2141.       lsrw    #4,d1
  2142.   |
  2143.   | Testing for NaN and Infinities
  2144.   |
  2145. !     movew    #0x7ff,d3
  2146. !     cmpw    d3,d0
  2147.       beq    0f
  2148. !     cmpw    d3,d1
  2149.       bne    nospec
  2150.       bra    1f
  2151.   |    first operand is special
  2152. ***************
  2153. *** 110,129 ****
  2154.       bne    retinf  |  Inf  * x == Inf
  2155.       orl    d5,d4
  2156.       beq    retnan  | Inf * 0 == NaN
  2157. ! retinf:    
  2158. !     clrl    d1
  2159. !     tstl    d2
  2160. !     bpl    0f
  2161. !     movel    #0xfff00000,d0
  2162.   return:    moveml    sp@+,d2-d7
  2163.       rts
  2164. ! 
  2165. ! 
  2166. ! 0:    movel    #0x7ff00000,d0
  2167. !     bra    return    
  2168. ! 
  2169. ! retnan: movel    #0x7fffffff,d0
  2170. !     moveql    #-1,d1
  2171.       bra    return
  2172.   |
  2173.   | v is special
  2174. --- 118,138 ----
  2175.       bne    retinf  |  Inf  * x == Inf
  2176.       orl    d5,d4
  2177.       beq    retnan  | Inf * 0 == NaN
  2178. ! |
  2179. ! |    Return Infinity with correct sign
  2180. ! |
  2181. ! retinf:    moveq    #0,d1
  2182. !     movel    #0xffe00000,d0    | we will return #0xfff00000 or #0x7ff00000
  2183. !     lslw    #1,d2
  2184. !     roxrl   #1,d0        | shift in high bit as given by d2
  2185.   return:    moveml    sp@+,d2-d7
  2186.       rts
  2187. ! |
  2188. ! |    Return NaN
  2189. ! |
  2190. ! retnan: moveql    #-1,d1
  2191. !     movel    d1,d0
  2192. !     bclr    #31,d0        | 0x7fffffff -> d0
  2193.       bra    return
  2194.   |
  2195.   | v is special
  2196. ***************
  2197. *** 137,159 ****
  2198.       bra    retinf
  2199.   |
  2200.   | end of NaN and Inf.
  2201. ! |    
  2202.   nospec:    subw    #16,sp        | multiplication accumulator
  2203.   
  2204. -     tstw    d0        | check for zero exponent - no leading "1"
  2205. -     beq    0f
  2206.       bset    #20,d6        | restore implied leading "1"
  2207. !     bra    1f
  2208. ! 0:    addw    #1,d0        | "normalize" exponent
  2209.   1:    movel    d6,d3
  2210.       orl    d7,d3
  2211.       beq    retz        | multiplying by zero
  2212.   
  2213. -     tstw    d1        | check for zero exponent - no leading "1"
  2214. -     beq    0f
  2215.       bset    #20,d4        | restore implied leading "1"
  2216. !     bra    1f
  2217. ! 0:    addw    #1,d1        | "normalize" exponent
  2218.   1:    movel    d4,d3
  2219.       orl    d5,d3
  2220.       beq    retz        | multiplying by zero
  2221. --- 146,168 ----
  2222.       bra    retinf
  2223.   |
  2224.   | end of NaN and Inf.
  2225. ! |
  2226.   nospec:    subw    #16,sp        | multiplication accumulator
  2227.   
  2228.       bset    #20,d6        | restore implied leading "1"
  2229. !     tstw    d0        | check for zero exponent - no leading "1"
  2230. !     bne    1f
  2231. !     bclr    #20,d6        | remove it
  2232. !     addw    #1,d0        | "normalize" exponent
  2233.   1:    movel    d6,d3
  2234.       orl    d7,d3
  2235.       beq    retz        | multiplying by zero
  2236.   
  2237.       bset    #20,d4        | restore implied leading "1"
  2238. !     tstw    d1        | check for zero exponent - no leading "1"
  2239. !     bne    1f
  2240. !     bclr    #20,d4        | remove it
  2241. !     addw    #1,d1        | "normalize" exponent
  2242.   1:    movel    d4,d3
  2243.       orl    d5,d3
  2244.       beq    retz        | multiplying by zero
  2245. ***************
  2246. *** 161,171 ****
  2247.       addw    d1,d0        | add exponents,
  2248.       subw    #BIAS8+16-11,d0    | remove excess bias, acnt for repositioning
  2249.   
  2250. !     clrl    sp@        | initialize 128-bit product to zero
  2251. !     clrl    sp@(4)
  2252. !     clrl    sp@(8)
  2253. !     clrl    sp@(12)
  2254. !     lea    sp@(8),a1    | accumulator pointer
  2255.   
  2256.   | see Knuth, Seminumerical Algorithms, section 4.3. algorithm M
  2257.   
  2258. --- 170,182 ----
  2259.       addw    d1,d0        | add exponents,
  2260.       subw    #BIAS8+16-11,d0    | remove excess bias, acnt for repositioning
  2261.   
  2262. !     lea    sp@,a1        | initialize 128-bit product to zero
  2263. !     moveq    #0,d3
  2264. !     movel    d3,a1@+
  2265. !     movel    d3,a1@+
  2266. !     movel    d3,a1@+
  2267. !     movel    d3,a1@
  2268. !     subqw    #4,a1        | an address of sp@(8) in a1
  2269.   
  2270.   | see Knuth, Seminumerical Algorithms, section 4.3. algorithm M
  2271.   
  2272. ***************
  2273. *** 202,209 ****
  2274.   
  2275.       moveml    sp@(2),d4-d7    | get the 112 valid bits
  2276.       clrw    d7        | (pad to 128)
  2277.   2:
  2278. !     cmpl    #0x0000ffff,d4    | multiply (shift) until
  2279.       bhi    3f        |  1 in upper 16 result bits
  2280.       cmpw    #9,d0        | give up for denormalized numbers
  2281.       ble    3f
  2282. --- 213,221 ----
  2283.   
  2284.       moveml    sp@(2),d4-d7    | get the 112 valid bits
  2285.       clrw    d7        | (pad to 128)
  2286. +     movel    #0x0000ffff,d3
  2287.   2:
  2288. !     cmpl    d3,d4        | multiply (shift) until
  2289.       bhi    3f        |  1 in upper 16 result bits
  2290.       cmpw    #9,d0        | give up for denormalized numbers
  2291.       ble    3f
  2292. ***************
  2293. *** 227,235 ****
  2294.       orb    #1,d1        | set "sticky bit" if any low-order set
  2295.   4:    addw    #16,sp        | remove accumulator from stack
  2296.       jmp    norm_df        | (result in d4/d5)
  2297. ! 
  2298. ! retz:    clrl    d0        | save zero as result
  2299. !     clrl    d1
  2300.       addw    #16,sp
  2301.       moveml    sp@+,d2-d7
  2302.       rts            | no normalizing neccessary
  2303. --- 239,249 ----
  2304.       orb    #1,d1        | set "sticky bit" if any low-order set
  2305.   4:    addw    #16,sp        | remove accumulator from stack
  2306.       jmp    norm_df        | (result in d4/d5)
  2307. ! |                | norm_df does not return here
  2308. ! retz:    moveq    #0,d0        | save zero as result
  2309. !     movel    d0,d1
  2310. !     lslw    #1,d2        | fill X bit
  2311. !     roxrl    #1,d0        | set high bit of d0 accordingly
  2312.       addw    #16,sp
  2313.       moveml    sp@+,d2-d7
  2314.       rts            | no normalizing neccessary
  2315. ===================================================================
  2316. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_mulsf3.cpp,v
  2317. retrieving revision 1.4
  2318. diff -c -r1.4 _mulsf3.cpp
  2319. *** 1.4    1993/03/01 17:34:14
  2320. --- _mulsf3.cpp    1993/05/29 20:28:36
  2321. ***************
  2322. *** 30,39 ****
  2323.       lea    0xfffffa50:w,a0
  2324.       movew    #0x4400,a0@(comm)    | load first argument to fp0
  2325.       cmpiw    #0x8900,a0@(resp)    | check
  2326. !     movel    a7@(4),a0@
  2327.       movew    #0x4427,a0@(comm)
  2328.       .long    0x0c688900, 0xfff067f8
  2329. !     movel    a7@(8),a0@
  2330.       movew    #0x6400,a0@(comm)    | result to d0
  2331.       .long    0x0c688900, 0xfff067f8
  2332.       movel    a0@,d0
  2333. --- 30,39 ----
  2334.       lea    0xfffffa50:w,a0
  2335.       movew    #0x4400,a0@(comm)    | load first argument to fp0
  2336.       cmpiw    #0x8900,a0@(resp)    | check
  2337. !     movel    sp@(4),a0@
  2338.       movew    #0x4427,a0@(comm)
  2339.       .long    0x0c688900, 0xfff067f8
  2340. !     movel    sp@(8),a0@
  2341.       movew    #0x6400,a0@(comm)    | result to d0
  2342.       .long    0x0c688900, 0xfff067f8
  2343.       movel    a0@,d0
  2344. ***************
  2345. *** 45,52 ****
  2346. --- 45,59 ----
  2347.   |
  2348.   | written by Kai-Uwe Bloem (I5110401@dbstu1.bitnet).
  2349.   | Based on a 80x86 floating point packet from comp.os.minix, written by P.Housel
  2350. + |
  2351. + | Revision 1.2.3 michal 05-93 (ntomczak@vm.ucs.ualberta.ca)
  2352. + |   + code smoothing
  2353. + |
  2354.   | patched by Olaf Flebbe (flebbe@tat.physik.uni-tuebingen.de)
  2355.   |
  2356. + | Revision 1.2.2 olaf 05-93
  2357. + |   + fixed a bug with -0.
  2358. + |
  2359.   | Revision 1.2.1 olaf 12-92:
  2360.   |   + added support for NaN and Infinites
  2361.   |   + added support for -0
  2362. ***************
  2363. *** 67,93 ****
  2364.       moveml    d2-d5,sp@-
  2365.       moveml    a0@,d4/d5    | d4 = v, d5 = u
  2366.   
  2367. ! 
  2368.       movel    d5,d0        | d0 = u.exp
  2369.       swap    d0
  2370.       movew    d0,d2        | d2 = u.sign
  2371.       lsrw    #7,d0
  2372. -     andw    #0xff,d0    | kill sign bit
  2373.   
  2374.       movel    d4,d1        | d1 = v.exp
  2375.       swap    d1
  2376. !     eorw    d1,d2        | d2 = u.sign ^ v.sign (in bit 31)
  2377.       lsrw    #7,d1
  2378. -     andw    #0xff,d1    | kill sign bit
  2379.   
  2380. -     andl    #0x7fffff,d5    | remove exponent from u.mantissa
  2381. -     andl    #0x7fffff,d4    | remove exponent from v.mantissa
  2382.   |
  2383.   | Testing for NaN and Infinities
  2384.   |
  2385. !     cmpb    #0xff,d0
  2386.       beq    0f
  2387. !     cmpb    #0xff,d1
  2388.       bne    nospec
  2389.       bra    1f
  2390.   |    first operand is special
  2391. --- 74,101 ----
  2392.       moveml    d2-d5,sp@-
  2393.       moveml    a0@,d4/d5    | d4 = v, d5 = u
  2394.   
  2395. !     movel    #0x7fffff,d3
  2396.       movel    d5,d0        | d0 = u.exp
  2397. +     andl    d3,d5        | remove exponent from u.mantissa
  2398.       swap    d0
  2399.       movew    d0,d2        | d2 = u.sign
  2400. +     bclr    #15,d0        | kill sign bit
  2401.       lsrw    #7,d0
  2402.   
  2403.       movel    d4,d1        | d1 = v.exp
  2404. +     andl    d3,d4        | remove exponent from v.mantissa
  2405.       swap    d1
  2406. !     eorw    d1,d2        | d2 = u.sign ^ v.sign (in bit 15)
  2407. !     bclr    #15,d1        | kill sign bit
  2408.       lsrw    #7,d1
  2409.   
  2410.   |
  2411.   | Testing for NaN and Infinities
  2412.   |
  2413. !     moveq    #-1,d3
  2414. !     cmpb    d3,d0
  2415.       beq    0f
  2416. !     cmpb    d3,d1
  2417.       bne    nospec
  2418.       bra    1f
  2419.   |    first operand is special
  2420. ***************
  2421. *** 100,116 ****
  2422.       tstl    d4
  2423.       beq    retnan  | Inf * 0 == NaN
  2424.   
  2425. ! retinf:    tstl    d2
  2426. !     bpl    0f
  2427. !     movel    #0xff800000,d0
  2428.   return:    moveml    sp@+,d2-d5
  2429.       rts
  2430. ! 
  2431. ! 
  2432. ! 0:    movel    #0x7f800000,d0
  2433. !     bra    return    
  2434. ! 
  2435. ! retnan: movel    #0x7fffffff,d0
  2436.       bra    return
  2437.   |
  2438.   | v is special
  2439. --- 108,126 ----
  2440.       tstl    d4
  2441.       beq    retnan  | Inf * 0 == NaN
  2442.   
  2443. ! |
  2444. ! |    Return Infinity with correct sign
  2445. ! |
  2446. ! retinf: movel    #0xff000000,d0  | we will return #0xff800000 or #0x7f800000
  2447. !     lslw    #1,d2
  2448. !     roxrl   #1,d0        | shift in high bit as given by d2
  2449.   return:    moveml    sp@+,d2-d5
  2450.       rts
  2451. ! |
  2452. ! |    Return NaN
  2453. ! |
  2454. ! retnan: moveql    #-1,d0
  2455. !     lsrl    #1,d0        | 0x7fffffff -> d0
  2456.       bra    return
  2457.   |
  2458.   | v is special
  2459. ***************
  2460. *** 124,144 ****
  2461.       bra    retinf
  2462.   |
  2463.   | end of NaN and Inf.
  2464. ! |    
  2465. ! nospec:    subw    #8,sp        | multiplication accumulator
  2466.       tstw    d0        | check for zero exponent - no leading "1"
  2467. !     beq    0f
  2468. !     bset    #23,d5        | restore implied leading "1"
  2469. !     bra    1f
  2470. ! 0:    addw    #1,d0        | "normalize" exponent
  2471.   1:    tstl    d5
  2472.       beq    retz        | multiplying zero
  2473.   
  2474. -     tstw    d1        | check for zero exponent - no leading "1"
  2475. -     beq    0f
  2476.       bset    #23,d4        | restore implied leading "1"
  2477. !     bra    1f
  2478. ! 0:    addw    #1,d1        | "normalize" exponent
  2479.   1:    tstl    d4
  2480.       beq    retz        | multiply by zero
  2481.   
  2482. --- 134,154 ----
  2483.       bra    retinf
  2484.   |
  2485.   | end of NaN and Inf.
  2486. ! |
  2487. ! nospec:    bset    #23,d5        | restore implied leading "1"
  2488. !     subw    #8,sp        | multiplication accumulator
  2489.       tstw    d0        | check for zero exponent - no leading "1"
  2490. !     bne    1f
  2491. !     bclr    #23,d5        | remove it
  2492. !     addqw    #1,d0        | "normalize" exponent
  2493.   1:    tstl    d5
  2494.       beq    retz        | multiplying zero
  2495.   
  2496.       bset    #23,d4        | restore implied leading "1"
  2497. !     tstw    d1        | check for zero exponent - no leading "1"
  2498. !     bne    1f
  2499. !     bclr    #23,d4        | remove it
  2500. !     addqw    #1,d1        | "normalize" exponent
  2501.   1:    tstl    d4
  2502.       beq    retz        | multiply by zero
  2503.   
  2504. ***************
  2505. *** 172,179 ****
  2506.                   | [TOP 16 BITS SHOULD BE ZERO !]
  2507.       moveml    sp@(2),d4-d5    | get the 48 valid mantissa bits
  2508.       clrw    d5        | (pad to 64)
  2509.   2:
  2510. !     cmpl    #0x0000ffff,d4    | multiply (shift) until
  2511.       bhi    3f        |  1 in upper 16 result bits
  2512.       cmpw    #9,d0        | give up for denormalized numbers
  2513.       ble    3f
  2514. --- 182,191 ----
  2515.                   | [TOP 16 BITS SHOULD BE ZERO !]
  2516.       moveml    sp@(2),d4-d5    | get the 48 valid mantissa bits
  2517.       clrw    d5        | (pad to 64)
  2518. + 
  2519. +     movel    #0x0000ffff,d3
  2520.   2:
  2521. !     cmpl    d3,d4        | multiply (shift) until
  2522.       bhi    3f        |  1 in upper 16 result bits
  2523.       cmpw    #9,d0        | give up for denormalized numbers
  2524.       ble    3f
  2525. ***************
  2526. *** 193,199 ****
  2527.   4:    addw    #8,sp        | remove accumulator from stack
  2528.       jmp    norm_sf        | (result in d4)
  2529.   
  2530. ! retz:    clrl    d0        | save zero as result
  2531.       addw    #8,sp
  2532.       moveml    sp@+,d2-d5
  2533.       rts            | no normalizing neccessary
  2534. --- 205,213 ----
  2535.   4:    addw    #8,sp        | remove accumulator from stack
  2536.       jmp    norm_sf        | (result in d4)
  2537.   
  2538. ! retz:    moveq    #0,d0        | save zero as result
  2539. !     lslw    #1,d2        | and set it sign as for d2
  2540. !     roxrl    #1,d0
  2541.       addw    #8,sp
  2542.       moveml    sp@+,d2-d5
  2543.       rts            | no normalizing neccessary
  2544. ===================================================================
  2545. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_normdf.cpp,v
  2546. retrieving revision 1.7
  2547. diff -c -r1.7 _normdf.cpp
  2548. *** 1.7    1992/12/28 07:56:08
  2549. --- _normdf.cpp    1993/05/29 20:28:37
  2550. ***************
  2551. *** 7,12 ****
  2552. --- 7,18 ----
  2553.   |
  2554.   | written by Kai-Uwe Bloem (I5110401@dbstu1.bitnet).
  2555.   | Based on a 80x86 floating point packet from comp.os.minix, written by P.Housel
  2556. + |
  2557. + | Revision 1.6.3 michal 05-93 (ntomczak@vm.ucs.ualberta.ca)
  2558. + |  + restored and ensured future synchronization with errno codes
  2559. + |  + removed bogus error when normalizing legitimate zero
  2560. + |  + small mods to shave off few cycles
  2561. + |
  2562.   | patched by Olaf Flebbe (flebbe@tat.physik.uni-tuebingen.de)
  2563.   |
  2564.   | Revision 1.6.2 olaf 12-92:
  2565. ***************
  2566. *** 73,84 ****
  2567.       movel    a0@+,d0        | get exponent
  2568.       movel    a0@+,d2        | get sign
  2569.       bpl    0f        | or bit 31 to bit 15 for later tests
  2570. !     orw    #0x8000,d2
  2571.   0:    movel    a0@+,d1        | rounding information
  2572.   
  2573. !     cmpl    #0x7fff,d0    | test exponent
  2574.       bgt    oflow
  2575. !     cmpl    #-0x8000,d0
  2576.       blt    retz
  2577.   # endif
  2578.   
  2579. --- 79,92 ----
  2580.       movel    a0@+,d0        | get exponent
  2581.       movel    a0@+,d2        | get sign
  2582.       bpl    0f        | or bit 31 to bit 15 for later tests
  2583. !     bset    #15,d2
  2584.   0:    movel    a0@+,d1        | rounding information
  2585.   
  2586. !     movel    #0x7fff,d3
  2587. !     cmpl    d3,d0        | test exponent
  2588.       bgt    oflow
  2589. !     notl    d3        | #-0x8000 -> d3
  2590. !     cmpl    d3,d0
  2591.       blt    retz
  2592.   # endif
  2593.   
  2594. ***************
  2595. *** 90,96 ****
  2596.       orl    d5,d3
  2597.       bne    1f
  2598.       tstb    d1
  2599. !     beq    retz
  2600.   1:
  2601.       movel    d4,d3
  2602.       andl    #0xfffff000,d3    | fast shift, 16 bits ?
  2603. --- 98,104 ----
  2604.       orl    d5,d3
  2605.       bne    1f
  2606.       tstb    d1
  2607. !     beq    retzok
  2608.   1:
  2609.       movel    d4,d3
  2610.       andl    #0xfffff000,d3    | fast shift, 16 bits ?
  2611. ***************
  2612. *** 122,128 ****
  2613.   4:
  2614.       andb    #1,d2
  2615.       orb    d2,d1        | make least sig bit sticky
  2616. !     movel    #0xfff00000,d6
  2617.   5:    movel    d4,d3        | multiply (shift) until
  2618.       andl    d6,d3        | one in implied position
  2619.       bne    6f
  2620. --- 130,136 ----
  2621.   4:
  2622.       andb    #1,d2
  2623.       orb    d2,d1        | make least sig bit sticky
  2624. !     asrl    #1,d6        | #0xfff00000 -> d6
  2625.   5:    movel    d4,d3        | multiply (shift) until
  2626.       andl    d6,d3        | one in implied position
  2627.       bne    6f
  2628. ***************
  2629. *** 137,143 ****
  2630.       bge    8f        | round down - no action neccessary
  2631.       negb    d1
  2632.       bvc    7f        | round up
  2633. !         movew   d5,d1           | tie case - round to even 
  2634.                                   | dont need rounding bits any more
  2635.           andw    #1,d1           | check if even
  2636.           beq     8f              | mantissa is even - no action necessary
  2637. --- 145,151 ----
  2638.       bge    8f        | round down - no action neccessary
  2639.       negb    d1
  2640.       bvc    7f        | round up
  2641. !         movew   d5,d1           | tie case - round to even
  2642.                                   | dont need rounding bits any more
  2643.           andw    #1,d1           | check if even
  2644.           beq     8f              | mantissa is even - no action necessary
  2645. ***************
  2646. *** 151,170 ****
  2647.       addw    #1,d0        | correct exponent for denormalized numbers
  2648.       bra    2b
  2649.   0:    movel    d4,d3        | check for rounding overflow
  2650. !     andl    #0xffe00000,d3
  2651.       bne    2b        | go back and renormalize
  2652.   8:
  2653.       movel    d4,d3        | check if normalization caused an underflow
  2654.       orl    d5,d3
  2655.       beq    retz
  2656. !     cmpw    #0,d0        | check for exponent overflow or underflow
  2657.       blt    retz
  2658.       cmpw    #2047,d0
  2659.       bge    oflow
  2660.   
  2661. !     lslw    #4,d0        | re-position exponent
  2662. !     andw    #0x8000,d2    | sign bit
  2663. !     orw    d2,d0
  2664.       swap    d0        | map to upper word
  2665.       clrw    d0
  2666.       andl    #0x0fffff,d4    | top mantissa bits
  2667. --- 159,179 ----
  2668.       addw    #1,d0        | correct exponent for denormalized numbers
  2669.       bra    2b
  2670.   0:    movel    d4,d3        | check for rounding overflow
  2671. !     asll    #1,d6        | #0xffe00000 -> d3
  2672. !     andl    d6,d3
  2673.       bne    2b        | go back and renormalize
  2674.   8:
  2675.       movel    d4,d3        | check if normalization caused an underflow
  2676.       orl    d5,d3
  2677.       beq    retz
  2678. !     tstw    d0        | check for exponent overflow or underflow
  2679.       blt    retz
  2680.       cmpw    #2047,d0
  2681.       bge    oflow
  2682.   
  2683. !     lslw    #5,d0        | re-position exponent - one bit too high
  2684. !     lslw    #1,d2        | get X bit
  2685. !     roxrw    #1,d0        | shift it into sign position
  2686.       swap    d0        | map to upper word
  2687.       clrw    d0
  2688.       andl    #0x0fffff,d4    | top mantissa bits
  2689. ***************
  2690. *** 175,187 ****
  2691.       rts
  2692.   
  2693.   retz:
  2694. !     moveq    #ERANGE,d0
  2695.       Emove   d0,Errno
  2696. !     clrl    d0        | return zero value
  2697. !     clrl    d1
  2698. !     tstw    d2
  2699. !     bpl    0f        | -0
  2700. !     bset    #31,d0
  2701.   0:    moveml    sp@+,d2-d7
  2702.       rts
  2703.   
  2704. --- 184,196 ----
  2705.       rts
  2706.   
  2707.   retz:
  2708. !     moveq    #Erange,d0
  2709.       Emove   d0,Errno
  2710. ! retzok:
  2711. !     moveq    #0,d0        | return zero value
  2712. !     movel    d0,d1
  2713. !     lslw    #1,d2        | set value of extension
  2714. !     roxrl    #1,d0        | and move it to hight bit of d0
  2715.   0:    moveml    sp@+,d2-d7
  2716.       rts
  2717.   
  2718. ***************
  2719. *** 191,206 ****
  2720.       pea    Stderr
  2721.       jbsr    _fprintf    |
  2722.       addql    #8,a7        |
  2723. !     moveq    #ERANGE,d0
  2724.       Emove    d0,Errno
  2725.   #endif    ERROR_CHECK
  2726.   
  2727.   |    moveml    pc@(__infinitydf),d0-d1 | return infinity value
  2728.       moveml    __infinitydf,d0-d1 | return infinty value
  2729. !     andw    #0x8000,d2    | get sign bit of argument
  2730. !     swap    d2
  2731. !     clrw    d2
  2732. !     orl    d2,d0
  2733.       moveml    sp@+,d2-d7    | should really cause trap ?!? (mjr: why?)
  2734.       rts
  2735.   
  2736. --- 200,215 ----
  2737.       pea    Stderr
  2738.       jbsr    _fprintf    |
  2739.       addql    #8,a7        |
  2740. !     moveq    #Erange,d0
  2741.       Emove    d0,Errno
  2742.   #endif    ERROR_CHECK
  2743.   
  2744.   |    moveml    pc@(__infinitydf),d0-d1 | return infinity value
  2745.       moveml    __infinitydf,d0-d1 | return infinty value
  2746. !     tstw    d2
  2747. !     bpl    1f
  2748. !     bset    #31,d0
  2749. ! 1:
  2750.       moveml    sp@+,d2-d7    | should really cause trap ?!? (mjr: why?)
  2751.       rts
  2752.   
  2753. ===================================================================
  2754. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_normsf.cpp,v
  2755. retrieving revision 1.5
  2756. diff -c -r1.5 _normsf.cpp
  2757. *** 1.5    1993/01/13 14:30:49
  2758. --- _normsf.cpp    1993/05/29 20:28:38
  2759. ***************
  2760. *** 7,12 ****
  2761. --- 7,18 ----
  2762.   |
  2763.   | written by Kai-Uwe Bloem (I5110401@dbstu1.bitnet).
  2764.   | Based on a 80x86 floating point packet from comp.os.minix, written by P.Housel
  2765. + |
  2766. + | Revision 1.4.4 michal 05-93 (ntomczak@vm.ucs.ualberta.ca)
  2767. + |  + restored and ensured future synchronization with errno codes
  2768. + |  + removed bogus error when normalizing legitimate zero
  2769. + |  + small mods to shave off few cycles
  2770. + |
  2771.   | patched by Olaf Flebbe (flebbe@tat.physik.uni-tuebingen.de)
  2772.   |
  2773.   | Revision 1.4.3 olaf 12-92:
  2774. ***************
  2775. *** 60,92 ****
  2776.   
  2777.       | internal entry for floating point package, saves time
  2778.       | d0=u.exp, d2=u.sign, d1=rounding bits, d4/d5=mantissa
  2779. !     | registers d2-d7 must be saved on the stack !
  2780.   norm_sf:
  2781.       tstl    d4        | rounding and u.mant == 0 ?
  2782.       bne    0f
  2783.       tstb    d1
  2784. !     beq    retz
  2785.   0:
  2786.       clrb    d2        | "sticky byte"
  2787. !     movel    #0xff000000,d5
  2788. ! 1:    tstw    d0        | divide (shift)
  2789.       ble    0f        |  denormalized number
  2790.       movel    d4,d3
  2791.       andl    d5,d3        |  or until no bits above 23
  2792.       beq    2f
  2793. ! 0:    addw    #1,d0        | increment exponent
  2794.       lsrl    #1,d4
  2795.       orb    d1,d2        | set "sticky"
  2796.       roxrb    #1,d1        | shift into rounding bits
  2797. !     bra    1b
  2798.   2:
  2799.       andb    #1,d2
  2800.       orb    d2,d1        | make least sig bit "sticky"
  2801. !     movel    #0xff800000,d5
  2802.   3:    movel    d4,d3        | multiply (shift) until
  2803.       andl    d5,d3        | one in "implied" position
  2804.       bne    4f
  2805. !     subw    #1,d0        | decrement exponent
  2806.       beq    4f        |  too small. store as denormalized number
  2807.       addb    d1,d1        | some doubt about this one *
  2808.       addxl    d4,d4
  2809. --- 66,98 ----
  2810.   
  2811.       | internal entry for floating point package, saves time
  2812.       | d0=u.exp, d2=u.sign, d1=rounding bits, d4/d5=mantissa
  2813. !     | registers d2-d5 must be saved on the stack !
  2814.   norm_sf:
  2815.       tstl    d4        | rounding and u.mant == 0 ?
  2816.       bne    0f
  2817.       tstb    d1
  2818. !     beq    retzok
  2819.   0:
  2820.       clrb    d2        | "sticky byte"
  2821. ! 1:    movel    #0xff000000,d5
  2822. ! 7:    tstw    d0        | divide (shift)
  2823.       ble    0f        |  denormalized number
  2824.       movel    d4,d3
  2825.       andl    d5,d3        |  or until no bits above 23
  2826.       beq    2f
  2827. ! 0:    addqw    #1,d0        | increment exponent
  2828.       lsrl    #1,d4
  2829.       orb    d1,d2        | set "sticky"
  2830.       roxrb    #1,d1        | shift into rounding bits
  2831. !     bra    7b
  2832.   2:
  2833.       andb    #1,d2
  2834.       orb    d2,d1        | make least sig bit "sticky"
  2835. !     asrl    #1,d5        | #0xff800000 -> d5
  2836.   3:    movel    d4,d3        | multiply (shift) until
  2837.       andl    d5,d3        | one in "implied" position
  2838.       bne    4f
  2839. !     subqw    #1,d0        | decrement exponent
  2840.       beq    4f        |  too small. store as denormalized number
  2841.       addb    d1,d1        | some doubt about this one *
  2842.       addxl    d4,d4
  2843. ***************
  2844. *** 98,104 ****
  2845.       bvc    5f        | round up
  2846.       movew   d4,d1           | tie case - round to even
  2847.                                   | dont need rounding bits any more
  2848. !     andw    #1,d1           | check if even        
  2849.       beq    6f              | mantissa is even - no action necessary
  2850.                                   | fall through
  2851.   
  2852. --- 104,110 ----
  2853.       bvc    5f        | round up
  2854.       movew   d4,d1           | tie case - round to even
  2855.                                   | dont need rounding bits any more
  2856. !     andw    #1,d1           | check if even
  2857.       beq    6f              | mantissa is even - no action necessary
  2858.                                   | fall through
  2859.   
  2860. ***************
  2861. *** 110,128 ****
  2862.       addw    #1,d0        | correct exponent for denormalized numbers
  2863.       bra    1b
  2864.   0:    movel    d4,d3        | check for rounding overflow
  2865. !     andl    #0xff000000,d3
  2866. !     bne    1b        | go back and renormalize
  2867.   6:
  2868.       tstl    d4        | check if normalization caused an underflow
  2869.       beq    retz
  2870. !     cmpw    #0,d0        | check for exponent overflow or underflow
  2871.       blt    retz
  2872.       cmpw    #255,d0
  2873.       bge    oflow
  2874.   
  2875. !     lslw    #7,d0        | re-position exponent
  2876. !     andw    #0x8000,d2    | sign bit
  2877. !     orw    d2,d0
  2878.       swap    d0        | map to upper word
  2879.       clrw    d0
  2880.       andl    #0x7fffff,d4    | top mantissa bits
  2881. --- 116,135 ----
  2882.       addw    #1,d0        | correct exponent for denormalized numbers
  2883.       bra    1b
  2884.   0:    movel    d4,d3        | check for rounding overflow
  2885. !     asll    #1,d5        | #0xff000000 -> d5
  2886. !     andl    d5,d3
  2887. !     bne    7b        | go back and renormalize
  2888.   6:
  2889.       tstl    d4        | check if normalization caused an underflow
  2890.       beq    retz
  2891. !     tstw    d0        | check for exponent overflow or underflow
  2892.       blt    retz
  2893.       cmpw    #255,d0
  2894.       bge    oflow
  2895.   
  2896. !     lslw    #8,d0        | re-position exponent - one bit too high
  2897. !     lslw    #1,d2        | get X bit
  2898. !     roxrw    #1,d0        | shift it into sign position
  2899.       swap    d0        | map to upper word
  2900.       clrw    d0
  2901.       andl    #0x7fffff,d4    | top mantissa bits
  2902. ***************
  2903. *** 130,163 ****
  2904.       moveml    sp@+,d2-d5
  2905.       rts
  2906.   
  2907. ! retz:    moveq    #ERANGE,d0
  2908.       Emove    d0,Errno
  2909. !     clrl    d0
  2910. !     tstw    d2
  2911. !     bpl    0f        | -0
  2912. !     bset    #31,d0
  2913. ! 0:    moveml    sp@+,d2-d5
  2914.       rts
  2915.   
  2916. ! oflow:    
  2917.   
  2918.   #ifdef    ERROR_CHECK
  2919. !     movel    d1,a7@-
  2920.       pea    pc@(LC0)
  2921.       pea    Stderr
  2922.       jbsr    _fprintf    |
  2923. !     addql    #8,a7        |
  2924. !     movel    a7@+,d1
  2925.   #endif    ERROR_CHECK
  2926.   
  2927.   |    movel    pc@(__infinitysf),d0    | return infinity value
  2928.       movel   __infinitysf,d0
  2929.   
  2930. !     btst    #15,d2        | transfer sign
  2931. !     beq    ofl_clear    | (mjr++)
  2932.       bset    #31,d0        |
  2933. -     moveml  sp@+,d2-d5
  2934. -     rts
  2935.   ofl_clear:
  2936.       moveml    sp@+,d2-d5    | should really cause trap ?!?
  2937.       rts
  2938. --- 137,168 ----
  2939.       moveml    sp@+,d2-d5
  2940.       rts
  2941.   
  2942. ! retz:    moveq    #Erange,d0
  2943.       Emove    d0,Errno
  2944. ! retzok:
  2945. !     moveq    #0,d0
  2946. !     lslw    #1,d2
  2947. !     roxrl    #1,d0        | sign of 0 is the same as of d2
  2948. !     moveml    sp@+,d2-d5
  2949.       rts
  2950.   
  2951. ! oflow:
  2952.   
  2953.   #ifdef    ERROR_CHECK
  2954. !     movel    d1,sp@-
  2955.       pea    pc@(LC0)
  2956.       pea    Stderr
  2957.       jbsr    _fprintf    |
  2958. !     addql    #8,sp        |
  2959. !     movel    sp@+,d1
  2960.   #endif    ERROR_CHECK
  2961.   
  2962.   |    movel    pc@(__infinitysf),d0    | return infinity value
  2963.       movel   __infinitysf,d0
  2964.   
  2965. !     tstw    d2        | transfer sign
  2966. !     bge    ofl_clear    | (mjr++)
  2967.       bset    #31,d0        |
  2968.   ofl_clear:
  2969.       moveml    sp@+,d2-d5    | should really cause trap ?!?
  2970.       rts
  2971. ===================================================================
  2972. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/_truncdf.cpp,v
  2973. retrieving revision 1.4
  2974. diff -c -r1.4 _truncdf.cpp
  2975. *** 1.4    1993/03/01 17:34:16
  2976. --- _truncdf.cpp    1993/05/29 20:28:40
  2977. ***************
  2978. *** 48,53 ****
  2979. --- 48,58 ----
  2980.   |
  2981.   | written by Kai-Uwe Bloem (I5110401@dbstu1.bitnet).
  2982.   | Based on a 80x86 floating point packet from comp.os.minix, written by P.Housel
  2983. + |
  2984. + | Revision 1.2.2 michal 05-93: (ntomczak@vm.ucs.ualberta.ca)
  2985. + |  + corrected bug in setting sign of returned Inf
  2986. + |  + code cleanup
  2987. + |
  2988.   | patched by Olaf Flebbe (flebbe@tat.physik.uni-tuebingen.de)
  2989.   |
  2990.   | Revision 1.2.1 olaf 12-92 :
  2991. ***************
  2992. *** 65,101 ****
  2993.   BIAS4    =    0x7F-1
  2994.   BIAS8    =    0x3FF-1
  2995.   
  2996. !     lea    sp@(4),a0    | parameter pointer
  2997.       moveml    d2-d5,sp@-    | save regs
  2998. -     moveml    a0@,d4-d5    | get number
  2999.   
  3000. !     movew    a0@,d0        | extract exponent
  3001.       movew    d0,d2        | extract sign
  3002.       lsrw    #4,d0
  3003. -     andw    #0x7ff,d0    | kill sign bit
  3004.       andl    #0x0fffff,d4    | remove exponent from mantissa
  3005.   |
  3006.   |
  3007.   |
  3008.       cmpw    #0x7ff,d0
  3009.       bne    nospec
  3010. !     orl    d4,d5
  3011.       bne    retnan
  3012. !     movel   #0x7f800000,d0
  3013. !     tstl    d4
  3014. !     bpl    return
  3015. !     bset    #31,d0
  3016.   return:    moveml    sp@+,d2-d5
  3017.       rts
  3018. ! retnan:    movel    #0x7fffffff,d0
  3019.       bra    return
  3020.   
  3021. ! nospec:    
  3022. !     tstw    d0        | check for zero exponent - no leading "1"
  3023. !     beq    0f        | for denormalized numbers
  3024.       bset    #20,d4        | restore implied leading "1"
  3025. !     bra    1f
  3026. ! 0:    addw    #1,d0        | "normalize" exponent
  3027.   1:
  3028.       addw    #BIAS4-BIAS8,d0    | adjust bias
  3029.   
  3030. --- 70,122 ----
  3031.   BIAS4    =    0x7F-1
  3032.   BIAS8    =    0x3FF-1
  3033.   
  3034. !     moveml    sp@(4),d0-d1    | get number
  3035.       moveml    d2-d5,sp@-    | save regs
  3036.   
  3037. !     movel    d0,d4        | save for norm_df
  3038. !     swap    d0        | extract exponent
  3039.       movew    d0,d2        | extract sign
  3040. +     bclr    #15,d0        | kill sign bit
  3041.       lsrw    #4,d0
  3042.       andl    #0x0fffff,d4    | remove exponent from mantissa
  3043.   |
  3044.   |
  3045.   |
  3046.       cmpw    #0x7ff,d0
  3047.       bne    nospec
  3048. !     orl    d4,d1
  3049.       bne    retnan
  3050. !     movel   #0xff00000,d0
  3051. !     lslw    #1,d2
  3052. !     roxrl    #1,d0
  3053.   return:    moveml    sp@+,d2-d5
  3054.       rts
  3055. ! retnan:    moveq    #-1,d0
  3056. !     lsrl    #1,d0        | #0x7fffffff -> d0
  3057.       bra    return
  3058.   
  3059. ! | Should we really return SNaN, which has a sign bit set accordingly??
  3060. ! | if yes, then the following code can be used instead -- mj
  3061. ! |
  3062. ! |    cmpw    #0x7ff,d0
  3063. ! |    bne    nospec
  3064. ! |    orl    d4,d1
  3065. ! |    bne    retnan
  3066. ! |    movel   #0xff00000,d0
  3067. ! |return:
  3068. ! |    lslw    #1,d2        | set X bit
  3069. ! |    roxrl    #1,d0        | roll in sign bit
  3070. ! |    moveml    sp@+,d2-d5
  3071. ! |    rts
  3072. ! |retnan:    moveq    #-1,d0
  3073. ! |    bra    return
  3074. ! 
  3075. ! nospec:    movel    d1,d5
  3076.       bset    #20,d4        | restore implied leading "1"
  3077. !     tstw    d0        | check for zero exponent - no leading "1"
  3078. !     bne    1f        | for denormalized numbers
  3079. !     bclr    #20,d4        | ... so do not do it but instead
  3080. !     addw    #1,d0        | "normalize" exponent
  3081.   1:
  3082.       addw    #BIAS4-BIAS8,d0    | adjust bias
  3083.   
  3084. ===================================================================
  3085. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/errbase.h,v
  3086. retrieving revision 1.1
  3087. diff -c -r1.1 errbase.h
  3088. *** 1.1    1992/10/09 20:35:29
  3089. --- errbase.h    1993/05/29 20:28:41
  3090. ***************
  3091. *** 2,12 ****
  3092.    * Errno defs to be included in *.ss floating point routines, from <errno.h>
  3093.    *   -- hyc@hanauma.jpl.nasa.gov, 9-17-92
  3094.    */
  3095.   
  3096.       .globl    _errno
  3097.   
  3098. ! EDOM    =    62
  3099. ! ERANGE    =    63
  3100.   
  3101.   #ifdef __MBASE__
  3102.   #define Errno    __MBASE__@(_errno)
  3103. --- 2,15 ----
  3104.    * Errno defs to be included in *.ss floating point routines, from <errno.h>
  3105.    *   -- hyc@hanauma.jpl.nasa.gov, 9-17-92
  3106.    */
  3107. + #define AssemB
  3108. + #include "errno.h"
  3109. + #undef  AssemB
  3110.   
  3111.       .globl    _errno
  3112.   
  3113. ! Edom    =    EDOM
  3114. ! Erange    =    ERANGE
  3115.   
  3116.   #ifdef __MBASE__
  3117.   #define Errno    __MBASE__@(_errno)
  3118. ===================================================================
  3119. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/frexp.cpp,v
  3120. retrieving revision 1.6
  3121. diff -c -r1.6 frexp.cpp
  3122. *** 1.6    1992/10/09 20:35:29
  3123. --- frexp.cpp    1993/05/29 20:28:43
  3124. ***************
  3125. *** 11,17 ****
  3126.    | ported to 68000 by Kai-Uwe Bloem, 12/89
  3127.    |  #1  original author: Peter S. Housel 9/21/88,01/17/89,03/19/89,5/24/89
  3128.    |  #2    added support for denormalized numbers            -kub-, 01/90
  3129. !  |  #   ported to gcc  ++jrb 04/90
  3130.    |-----------------------------------------------------------------------------
  3131.   
  3132.   BIAS8    =    0x3ff - 1
  3133. --- 11,19 ----
  3134.    | ported to 68000 by Kai-Uwe Bloem, 12/89
  3135.    |  #1  original author: Peter S. Housel 9/21/88,01/17/89,03/19/89,5/24/89
  3136.    |  #2    added support for denormalized numbers            -kub-, 01/90
  3137. !  |  #3  ported to gcc  ++jrb 04/90
  3138. !  |  #4  support for signed 0 - consistently with the rest - michal, 05/93
  3139. !  |      and some code cleanup
  3140.    |-----------------------------------------------------------------------------
  3141.   
  3142.   BIAS8    =    0x3ff - 1
  3143. ***************
  3144. *** 25,62 ****
  3145.   #else
  3146.       clrl    a0@
  3147.   #endif
  3148. !     tstl    sp@(4)        | 1st arg == 0 ?
  3149.       bne    nonzero
  3150. !     tstl    sp@(8)
  3151.       bne    nonzero
  3152. !     clrl    d0
  3153. !     clrl    d1
  3154. !     rts
  3155.   nonzero:
  3156. !     lea    sp@(4),a1
  3157. !     moveml    d2-d7,sp@-
  3158. ! 2:
  3159. !     movew    a1@,d0        | extract value.exp
  3160. !     movew    d0,d2        | extract value.sign    ++jrb
  3161.       lsrw    #4,d0
  3162. !     andw    #0x7ff,d0    | kill sign bit
  3163.       cmpw    #BIAS8,d0    | get out of loop if finally (a1) in [0.5,1.0)
  3164.       beq    3f
  3165.   
  3166. !     andw    #0x0f,a1@    | remove exponent from value.mantissa
  3167. !     tstw    d0        | check for zero exponent - no leading 1
  3168. !     beq    0f
  3169. !     orw    #0x10,a1@    | restore implied leading 1
  3170. !     bra    1f
  3171. ! 0:    addw    #1,d0
  3172. ! 1:
  3173.       movel    a1@,d1        | check for zero
  3174.       orl    a1@(4),d1
  3175.       beq    3f        | if zero, all done : exp = 0, num = 0.0
  3176. ! 
  3177.       subw    #BIAS8,d0    | remove bias
  3178. ! #ifdef __MSHORT__        
  3179. !     addw    d0,a0@        | add current exponent in
  3180.   #else
  3181.       extl    d0
  3182.       addl    d0,a0@        | add current exponent in
  3183. --- 27,69 ----
  3184.   #else
  3185.       clrl    a0@
  3186.   #endif
  3187. !     lea    sp@(4),a1    | sp@(4) -> a1 - an address of value
  3188. ! |                |   and a1@(8) - an address of *eptr
  3189. !     movel    a1@,d0        | keep value.sign
  3190. !     movel    d0,d1
  3191. !     bclr    #31,d1        | kill sign bit
  3192. !     tstl    d1        | 1st arg == 0 ?
  3193.       bne    nonzero
  3194. !     tstl    a1@(4)
  3195.       bne    nonzero
  3196. !     rts            | if 0 or -0 then d0 and d1 already set
  3197.   nonzero:
  3198. !     movel    d2,sp@-
  3199. ! 2:                | return here when looping
  3200. !     swap    d0        | sign and exponent into lower 16 bits
  3201. !     movew    d0,d2        | set d2 for norm_df
  3202. !     bclr    #15,d0        | kill sign bit
  3203.       lsrw    #4,d0
  3204. ! 
  3205.       cmpw    #BIAS8,d0    | get out of loop if finally (a1) in [0.5,1.0)
  3206.       beq    3f
  3207.   
  3208. !     moveq    #0x0f,d1    | remove exponent from value.mantissa
  3209. !     andb    d2,d1        | four upper bits of value in d1
  3210. !     bset    #4,d1        | implied leading 1
  3211. !     tstw    d0        | check for zero exponent
  3212. !     bne    1f
  3213. !     addqw    #1,d0
  3214. !     bclr    #4,d1        | nah, we do not need stinkin leadin 1
  3215. ! 1:    movew    d1,a1@        | save results of our efforts
  3216. ! 
  3217.       movel    a1@,d1        | check for zero
  3218.       orl    a1@(4),d1
  3219.       beq    3f        | if zero, all done : exp = 0, num = 0.0
  3220. !                 | sign of zero is correct
  3221.       subw    #BIAS8,d0    | remove bias
  3222. ! #ifdef __MSHORT__
  3223. !     addw    d0,a0        | add current exponent in
  3224.   #else
  3225.       extl    d0
  3226.       addl    d0,a0@        | add current exponent in
  3227. ***************
  3228. *** 68,79 ****
  3229.       moveml    d2-d7,sp@-    | ... need to copy with -mshort)
  3230.       moveml    a1@,d4-d5
  3231.       jmp    norm_df        | normalize result
  3232. ! L0:
  3233.       moveml    d0-d1,a1@
  3234.       bra    2b        | loop around to catch denormalized numbers
  3235.   3:
  3236.       moveml    a1@,d0-d1
  3237. !     moveml    sp@+,d2-d7
  3238.                   | d0-d1 has normalized mantissa
  3239.       rts
  3240.   
  3241. --- 75,87 ----
  3242.       moveml    d2-d7,sp@-    | ... need to copy with -mshort)
  3243.       moveml    a1@,d4-d5
  3244.       jmp    norm_df        | normalize result
  3245. ! L0:                | norm_df does not affect a0 or a1
  3246. ! |                | but it pops d2-d7
  3247.       moveml    d0-d1,a1@
  3248.       bra    2b        | loop around to catch denormalized numbers
  3249.   3:
  3250.       moveml    a1@,d0-d1
  3251. !     movel    sp@+,d2
  3252.                   | d0-d1 has normalized mantissa
  3253.       rts
  3254.   
  3255. ***************
  3256. *** 94,100 ****
  3257.   |#######################################################################
  3258.   | a0:    pointer to n
  3259.   |     sp@(12) address of resulting exponent (n)
  3260. !     
  3261.       .text; .even
  3262.   
  3263.   .globl _frexp
  3264. --- 102,108 ----
  3265.   |#######################################################################
  3266.   | a0:    pointer to n
  3267.   |     sp@(12) address of resulting exponent (n)
  3268. ! 
  3269.       .text; .even
  3270.   
  3271.   .globl _frexp
  3272. ***************
  3273. *** 161,171 ****
  3274.   |    bmi.b    www
  3275.   | is coded by
  3276.   |    .word    0x4a68,0xfff0,0x6bfa        | test
  3277. !     
  3278.       .text; .even
  3279.   
  3280.   .globl _frexp
  3281. !     
  3282.       tstl    sp@(4)        | 1st arg == 0 ?
  3283.       bne    nonzero
  3284.       tstl    sp@(8)
  3285. --- 169,179 ----
  3286.   |    bmi.b    www
  3287.   | is coded by
  3288.   |    .word    0x4a68,0xfff0,0x6bfa        | test
  3289. ! 
  3290.       .text; .even
  3291.   
  3292.   .globl _frexp
  3293. ! 
  3294.       tstl    sp@(4)        | 1st arg == 0 ?
  3295.       bne    nonzero
  3296.       tstl    sp@(8)
  3297. ===================================================================
  3298. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/ldexp.cpp,v
  3299. retrieving revision 1.6
  3300. diff -c -r1.6 ldexp.cpp
  3301. *** 1.6    1992/10/09 20:35:29
  3302. --- ldexp.cpp    1993/05/29 20:28:46
  3303. ***************
  3304. *** 42,47 ****
  3305. --- 42,48 ----
  3306.    |  #2    added support for denormalized numbers            -kub-, 01/90
  3307.    |  #3  ported to gcc  ++jrb  04/90
  3308.    |  #4  handle exponent overflow when ints are 32 bits        -kub-, 04/90
  3309. +  |  #5  make returned zero signed                  mj,  05/93
  3310.    |-----------------------------------------------------------------------------
  3311.   
  3312.       lea    sp@(4),a1
  3313. ***************
  3314. *** 49,64 ****
  3315.   
  3316.       movew    a1@,d0        | extract value.exp
  3317.       movew    d0,d2        | extract value.sign
  3318.       lsrw    #4,d0
  3319. -     andw    #0x7ff,d0    | kill sign bit
  3320.   
  3321. !     andw    #0x0f,a1@    | remove exponent from value.mantissa
  3322. !     tstw    d0        | check for zero exponent - no leading 1
  3323. !     beq    0f
  3324. !     orw    #0x10,a1@    | restore implied leading 1
  3325. !     bra    1f
  3326. ! 0:    addw    #1,d0
  3327. ! 1:
  3328.   #ifdef __MSHORT__
  3329.       addw    a1@(8),d0    | add in exponent
  3330.       extl    d0
  3331. --- 50,66 ----
  3332.   
  3333.       movew    a1@,d0        | extract value.exp
  3334.       movew    d0,d2        | extract value.sign
  3335. +     bclr    #15,d0        | kill sign bit
  3336.       lsrw    #4,d0
  3337.   
  3338. !     moveq    #0x0f,d3    | remove exponent from value.mantissa
  3339. !     andb    d2,d3        | four upper bits of value in d3
  3340. !     bset    #4,d3        | implied leading 1
  3341. !     tstw    d0        | check for zero exponent
  3342. !     bne    1f
  3343. !     addqw    #1,d0
  3344. !     bclr    #4,d3        | nah, we do not need stinkin leadin 1
  3345. ! 1:    movew    d3,a1@        | save results of our efforts
  3346.   #ifdef __MSHORT__
  3347.       addw    a1@(8),d0    | add in exponent
  3348.       extl    d0
  3349. ***************
  3350. *** 68,74 ****
  3351.   #endif
  3352.       cmpl    #-53,d0        | hmm. works only if 1 in implied position...
  3353.       ble    retz        | range error - underflow
  3354. !     cmpl    #2047,d0
  3355.       bge    rangerr        | range error - overflow
  3356.   
  3357.       clrw    d1        | zero rounding bits
  3358. --- 70,76 ----
  3359.   #endif
  3360.       cmpl    #-53,d0        | hmm. works only if 1 in implied position...
  3361.       ble    retz        | range error - underflow
  3362. !     cmpl    #0x7ff,d0
  3363.       bge    rangerr        | range error - overflow
  3364.   
  3365.       clrw    d1        | zero rounding bits
  3366. ***************
  3367. *** 77,89 ****
  3368.   
  3369.   retz:
  3370.       moveq    #0,d0        | zero return value
  3371. !     moveq    #0,d1
  3372.       jra    L0
  3373.   
  3374.   rangerr:
  3375.   
  3376.   #ifdef ERROR_CHECK
  3377. !     moveq    #ERANGE,d0
  3378.       Emove    d0,Errno
  3379.       pea    pc@(_Overflow)    | for printf
  3380.       pea    Stderr        |
  3381. --- 79,93 ----
  3382.   
  3383.   retz:
  3384.       moveq    #0,d0        | zero return value
  3385. !     movel    d0,d1
  3386. !     lslw    #1,d2        | transfer argument sign
  3387. !     roxrl    #1,d0
  3388.       jra    L0
  3389.   
  3390.   rangerr:
  3391.   
  3392.   #ifdef ERROR_CHECK
  3393. !     moveq    #Erange,d0
  3394.       Emove    d0,Errno
  3395.       pea    pc@(_Overflow)    | for printf
  3396.       pea    Stderr        |
  3397. ***************
  3398. *** 92,101 ****
  3399.   #endif ERROR_CHECK
  3400.   
  3401.       moveml    __infinitydf,d0-d1 | return HUGE_VAL (same as in <math.h>)
  3402. !     andw    #0x8000,d2    | get sign bit of argument
  3403. !     swap    d2
  3404. !     clrw    d2
  3405. !     orl    d2,d0
  3406.   L0:
  3407.       moveml    sp@+,d2-d7    | pop saved reggies
  3408.       rts
  3409. --- 96,104 ----
  3410.   #endif ERROR_CHECK
  3411.   
  3412.       moveml    __infinitydf,d0-d1 | return HUGE_VAL (same as in <math.h>)
  3413. !     tstw    d2
  3414. !     bge    L0
  3415. !     bset    #31,d0
  3416.   L0:
  3417.       moveml    sp@+,d2-d7    | pop saved reggies
  3418.       rts
  3419. ***************
  3420. *** 201,212 ****
  3421.   error_plus:
  3422.       swap    d0
  3423.       moveml    d0-d1,a7@-    | print error message
  3424. !     moveq    #ERANGE,d0    | Overflow: errno = ERANGE
  3425.       Emove    d0,Errno
  3426.       pea    pc@(_Domain)    | for printf
  3427.       bra    error_exit    |
  3428.   error_nan:
  3429. !     moveq    #EDOM,d0    | NAN => errno = EDOM
  3430.       Emove    d0,Errno
  3431.       moveml    a0@(24),d0-d1    | result = +inf
  3432.       moveml    d0-d1,a7@-    | print error message
  3433. --- 204,215 ----
  3434.   error_plus:
  3435.       swap    d0
  3436.       moveml    d0-d1,a7@-    | print error message
  3437. !     moveq    #Erange,d0    | Overflow: errno = ERANGE
  3438.       Emove    d0,Errno
  3439.       pea    pc@(_Domain)    | for printf
  3440.       bra    error_exit    |
  3441.   error_nan:
  3442. !     moveq    #Edom,d0    | NAN => errno = EDOM
  3443.       Emove    d0,Errno
  3444.       moveml    a0@(24),d0-d1    | result = +inf
  3445.       moveml    d0-d1,a7@-    | print error message
  3446. ===================================================================
  3447. RCS file: /net/acae127/home/bammi/etc/src/master/atari/co/modf.cpp,v
  3448. retrieving revision 1.2
  3449. diff -c -r1.2 modf.cpp
  3450. *** 1.2    1992/10/09 20:35:29
  3451. --- modf.cpp    1993/05/29 20:28:48
  3452. ***************
  3453. *** 1,18 ****
  3454.   #if !defined (__M68881__) && !defined (sfp004)
  3455.   
  3456.    | take floating point to integer and fractional pieces
  3457. !  | 
  3458.    | C interface
  3459.    |  double modf( double value, double *iptr)
  3460.    |  returns fractional part of value
  3461.    |       in *iptr returns the integral part
  3462.    |       such that (*iptr + fractional) == value
  3463. !  | 
  3464.    |-----------------------------------------------------------------------------
  3465.    | ported to 68000 by Kai-Uwe Bloem, 12/89
  3466.    |  #1  original author: Peter S. Housel 9/21/88,01/17/89,03/19/89,5/24/89
  3467.    |  #2  replaced shifts by swap if possible for speed increase    -kub-, 01/90
  3468.    |  #3  ported to gcc ++jrb 03/90
  3469.    |-----------------------------------------------------------------------------
  3470.   
  3471.   BIAS8    =    0x3ff - 1
  3472. --- 1,20 ----
  3473.   #if !defined (__M68881__) && !defined (sfp004)
  3474.   
  3475.    | take floating point to integer and fractional pieces
  3476. !  |
  3477.    | C interface
  3478.    |  double modf( double value, double *iptr)
  3479.    |  returns fractional part of value
  3480.    |       in *iptr returns the integral part
  3481.    |       such that (*iptr + fractional) == value
  3482. !  |
  3483.    |-----------------------------------------------------------------------------
  3484.    | ported to 68000 by Kai-Uwe Bloem, 12/89
  3485.    |  #1  original author: Peter S. Housel 9/21/88,01/17/89,03/19/89,5/24/89
  3486.    |  #2  replaced shifts by swap if possible for speed increase    -kub-, 01/90
  3487.    |  #3  ported to gcc ++jrb 03/90
  3488. +  |  #4  replaced by a completely new, smaller and faster implementation,
  3489. +  |        Michal Jaegermann, ntomczak@vm.ucs.ualberta.ca - 05/93
  3490.    |-----------------------------------------------------------------------------
  3491.   
  3492.   BIAS8    =    0x3ff - 1
  3493. ***************
  3494. *** 21,104 ****
  3495.       .globl _modf
  3496.   _modf:
  3497.       lea    sp@(4),a0    | a0 -> float argument
  3498. -     movel    sp@(12),a1    | a1 -> ipart result
  3499.       moveml    d2-d7,sp@-    | save d2-d7
  3500.   
  3501. !     movew    a0@,d0        | extract value.exp
  3502. !     movew    d0,d2        | extract value.sign
  3503. !     lsrw    #4,d0
  3504. !     andw    #0x7ff,d0    | kill sign bit
  3505. ! 
  3506. !     cmpw    #BIAS8,d0
  3507. !     bge    1f        | fabs(value) >= 1.0
  3508. ! 
  3509. !     clrl    a1@        | store zero as the integer part
  3510. !     clrl    a1@(4)
  3511. !     moveml    a0@,d0-d1    | return entire value as fractional part
  3512. !     jra    L0
  3513.   
  3514.   1:
  3515. !     cmpw    #BIAS8+53,d0    | all integer, with no fractional part ?
  3516. !     blt    2f        | no, mixed
  3517.   
  3518. !     movel    a0@,a1@        | store entire value as the integer part
  3519. !     movel    a0@(4),a1@(4)
  3520.       moveq    #0,d0        | return zero as fractional part
  3521. !     moveq    #0,d1
  3522. ! L0:
  3523. !     moveml    sp@+,d2-d7    | restore saved d2-d7
  3524. !     rts
  3525.   2:
  3526. ! |    moveml    d3-d7,sp@-    | save some registers (d2 is pushed already)
  3527. !     moveml    a0@,d4-d5    | get value
  3528. ! 
  3529. !     andl    #0x0fffff,d4    | remove exponent from value.mantissa
  3530. !     orl    #0x100000,d4    | restore implied leading "1"
  3531.   
  3532. !     moveq    #0,d6        | zero fractional part
  3533. !     moveq    #0,d7
  3534.   3:
  3535. !     cmpw    #BIAS8+37,d0    | fast shift, 16 bits ?
  3536. !     bgt    4f
  3537. !     movew    d6,d7        | shift down 16 bits
  3538. !     movew    d5,d6
  3539. !     movew    d4,d5
  3540. !     clrw    d4
  3541. !     swap    d7
  3542. !     swap    d6
  3543. !     swap    d5
  3544. !     swap    d4
  3545. !     addw    #16,d0
  3546. !     bra    3b
  3547.   4:
  3548. !     cmpw    #BIAS8+53,d0    | done ?
  3549. !     bge    5f
  3550. !     lsrl    #1,d4        | shift integer part
  3551. !     roxrl    #1,d5
  3552. ! 
  3553. !     roxrl    #1,d6        | shift high bit into fractional part
  3554. !     roxrl    #1,d7
  3555. ! 
  3556. !     addw    #1,d0        | increment ipart exponent
  3557. !     bra    4b        | keep shifting
  3558. ! 5:
  3559. !     | normalize ipart (all values are in correct reggies)
  3560. !      | save  a1, save d2-d7 that norm_df will pop
  3561. !     movel    a1,sp@-
  3562. !     pea    L1        | set up return address
  3563. !     moveml    d2-d7,sp@-    | norm_df will pop this
  3564. !     clrw    d1
  3565. !     jmp    norm_df        | go do it
  3566. ! L1:                | norm_df will rts to here
  3567. !     movel    sp@+,a1        | pop saved a1
  3568. !     moveml    d0-d1,a1@    | store result into ipart
  3569. !     
  3570. !     | norm fractional part
  3571. !     movel    d6,d4        | get frac into d4/d5
  3572. !     movel    d7,d5
  3573.       clrw    d1        | rounding = 0
  3574. !     movew    #BIAS8-11,d0    | set frac part exponent, sign already in d2
  3575. !     jmp    norm_df        | norm_df will pop d2/d7 we save before
  3576.                   | it will return to our caller via rts
  3577.                   | with result in d0-d1
  3578.   
  3579. --- 23,88 ----
  3580.       .globl _modf
  3581.   _modf:
  3582.       lea    sp@(4),a0    | a0 -> float argument
  3583.       moveml    d2-d7,sp@-    | save d2-d7
  3584. +     moveml    a0@+,d0-d1
  3585. +     movel   a0@,a1        | a1 -> ipart result
  3586.   
  3587. !     movel    d0,d2        | calculate exponent
  3588. !     swap    d2
  3589. !     bclr    #15,d2        | kill sign bit
  3590. !     lsrw    #4,d2        | exponent in lower 12 bits of d2
  3591. ! 
  3592. !     cmpw    #BIAS8,d2
  3593. !     bgt    1f        | fabs(value) >= 1.0
  3594. ! |                | return entire value as fractional part
  3595. !     clrl    a1@+        | d0, d1 already ok
  3596. !     clrl    a1@        | make integer part 0
  3597.   
  3598. + 0:
  3599. +     moveml    sp@+,d2-d7    | restore saved d2-d7
  3600. +     rts
  3601. + 
  3602.   1:
  3603. !     movew    #BIAS8+53,d3
  3604. !     subw    d2,d3        | compute position of "binary point"
  3605. !     bgt    2f        | branch if we do have fractional part
  3606.   
  3607. !     moveml  d0-d1,a1@    | store entire value as the integer part
  3608.       moveq    #0,d0        | return zero as fractional part
  3609. !     movel    d0,d1
  3610. !     jra    0b
  3611.   2:
  3612. !     movel    d1,d5        | save for computation of fractional part
  3613.   
  3614. !     moveq    #32,d6
  3615. !     cmpw    d6,d3
  3616. !     blt    3f        | jump if "binary point" in a lower part
  3617. !     movel    d0,d4
  3618. !     subw    d6,d3
  3619. !     moveq    #0,d6        | compute mask for splitting
  3620. !     bset    d3,d6
  3621. !     negl    d6
  3622. !     andl    d6,d0        | this is integer part
  3623. !     moveq    #0,d1
  3624. !     notl    d6
  3625. !     andl    d6,d4        | and denormalized fractional part
  3626. !     jra    4f
  3627.   3:
  3628. !     moveq    #0,d6        | splitting on lower part
  3629. !     bset    d3,d6
  3630. !     negl    d6
  3631. !     andl    d6,d1        | this is integer part
  3632. !     moveq    #0,d4        | nothing in an upper fraction
  3633. !     notl    d6
  3634. !     andl    d6,d5        | and clear those unneded bits
  3635.   4:
  3636. !     moveml    d0-d1,a1@    | store computed integer part
  3637. ! 
  3638. !     swap    d0
  3639. !     exg    d0,d2        | set registers for norm_df
  3640.       clrw    d1        | rounding = 0
  3641. ! |                | normalize fractional part
  3642. !     jmp    norm_df        | norm_df will pop d2/d7 we saved before
  3643.                   | it will return to our caller via rts
  3644.                   | with result in d0-d1
  3645.   
  3646. ***************
  3647. *** 109,115 ****
  3648.   |       manually optimized by Michael Ritzert
  3649.   |
  3650.   |    double modf( double X, double * IP )
  3651. ! |    
  3652.   |    30.11.92
  3653.   |    ritzert@dfg.dbp.de
  3654.   |
  3655. --- 93,99 ----
  3656.   |       manually optimized by Michael Ritzert
  3657.   |
  3658.   |    double modf( double X, double * IP )
  3659. ! |
  3660.   |    30.11.92
  3661.   |    ritzert@dfg.dbp.de
  3662.   |
  3663.